Exact solution of the Falicov-Kimball model with dynamical mean-field theory 
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The Falicov-Kimball model was introduced in 1969 as a statistical model for metal-insulator tran- 
sitions; it includes itinerant and localized electrons that mutually interact with a local Coulomb 
interaction and is the simplest model of electron correlations. It can be solved exactly with dy- 
namical mean-field theory in the limit of large spatial dimensions which provides an interesting 
benchmark for the physics of locally correlated systems. In this review, we develop the formalism 
for solving the Falicov-Kimball model from a path-integral perspective, and provide a number of 
expressions for single and two-particle properties. We examine many important theoretical results 
that show the absence of fermi-liquid features and provide a detailed description of the static and 
dynamic correlation functions and of transport properties. The parameter space is rich and one 
finds a variety of many-body features like metal-insulator transitions, classical valence fluctuating 
transitions, metamagnetic transitions, charge density wave order-disorder transitions, and phase 
separation. At the same time, a number of experimental systems have been discovered that show 
anomalies related to Falicov-Kimball physics [including YbInCu4, EuNi2(Sii— !Gei)2, N1I2 and 
Ta^N]. 



PACS numbers: 71.10.-W, 71.30.+h, 71.45.-d, 72.10.-d 
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A. Brief History 



The Falicov-Kimball model ijFalicov and Kimball 

1969) was introduced in 1969 to describe metal-insulator 
transitions in a number of rare-earth and transi tion-metal 
comp ounds [but see Hubbard's earlier work ijHubbardl 
Il963|) where the spinless version of the Falicov-Kimball 
model was introduced as an approximate solution to 
the Hubbard model; one assumes that one species 
of spin does not hop and is frozen on the lattice]. 
The initial work by Falicov and collaborators focused 
primarily on analyzing the thermodynamics of the 
metal-insulat or transition with a static mean-field the- 
ory a pproach l(Falicov and Kimbalill96lilR,amirez et all 

1970) . The resulting solutions displayed both continu- 
ous and discontinuous metal-insulator phase transitions, 
and they could fit the conductivity of a wide variety 
of transition-metal and rare-earth compounds with their 
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results. Next, iRamirez and Falicovl l)l97ll) applied the 
model to describe the a — 7 phase transition in Cerium. 
Again, a number of thermodynamic quantities were ap- 
proximated well by the model, but it did not examine any 
effects associated with Kondo screening of the /-electrons 
[and subsequently was discarded in favor of the Kondo 
volume collapse picture ijAllen and Martini Il9 82)]. 

Interest in the model waned once iPlischkel lll972l) 
showed that when the coherent-potentia l-approximation 
(CPA) l|Soverl Il967t IVelickv et a ll Il968h was applied to 
it, all first-order phase transitions disappeared and the 
solutions only displayed smooth crossovers from a metal 
to an insulat or [this claim was strongly refuted by Fal- 
icov's group ijGoncalves da Silva and Falicovl Il972|) but 
interest in the model was limited for almost 15 years]. 

The field was revitalized by mathematical physicists in 
the mid 1980's, who realized that the spinless version of 
this model is the simplest correlated system that displays 
long-range order at low temperatures and for dimensions 
greater than one. Indeed, two gro ups produced indepen- 
dent proofs of the long-range order {B randt a nd Schmidt! 
11 9Sd 11 9871 l kennedv and L iehi ll 9Sfit iLiebl H 9Sfit> . The 
work by Kennedy and Lieb rediscovered Hubbard's 
original approximation that yields the spinless Falicov- 
Kimball model, and also provided a new interpretation 
of the model for the physics of crystallization. A number 
of other exact results followed including (i) a proof of 
no quantum-mechanical mi xed valence (or spontaneous 
hybri dization) at finite T ijSubrahmanvam and Barmal 
1988) [based on the prese nce of a local gauge symme- 
try and Elitzur's theorem ijElitzurl Il975|) ]: (ii) proofs of 
phase separation and of periodic orde ring in one dimen - 
sion (with large interaction strength) ijLembergerl tl992) : 
(iii) proofs about ground-state properties i n two dimen- 
sions (also at large i nteraction strength) llHallerl l2000l 
lHaller and Kennedvl l200ll iKennedvl 11994 Il998|k (iv) a 
proof of phase separ ation in one dimensio n and small 
interaction strength l|Freericks et all Il996j) : and (v) a 
proof of phase sepa ration for large i nteractio n strength 
and all dimensions l)Freericks et all l2002albD . Most of 
these rig orous re s ults have already been summariz ed in 
reviews l)Gruberl Il999l iGruber and Macrisl Il996|) . In 
addition, a series of numerical c alculations were per- 
formed in one and two dimensions llde Vries et a H 119931 



forme d m one and two dimensions Ide^Hes^^zL 
1994 IFreericks and Falicovl Il99fi IGruber et all 



1994: 



Michielsenl Il993t IWatson and Lemanskil Il995[) While 



not providing complete results for the model, the nu- 
merics do illustrate a number of important trends in the 
physics of the FK model. 

At about the same time, there was a parallel devel- 
opment of the dynamical mean-field theory (DMFT), 
which is what we concentrat e on in this review. 
The D MFT was invented by iMetzner and Vollhardtl 
(1989). Almost immediately after the idea that in 
large spatial dimensions the self energy becomes lo- 
cal, Brandt and collaborators showed how to solve 
the static problem exactly (requiring no quantum 
Monte Carlo), thereby providing the exact solution of 



the Falicov-Kimball model (Brand^and Flcddcriohann, 
1992 f IBrandt et all [l990t Brandt and Mielschl Il98£ , 
19901 119911 IBrandt and Urbanekl Il992|) . This work is 



the extension of Onsager's famous solution for the tran- 
sition temperature of the two-dimensional Ising model to 
the fermionic case (and large dimensions). These series of 
papers revolutionized Falicov-Kimball model physics and 
provided the only exact quantitative results for electronic 
phase transitions in the thermodynamic limit for all val- 
ues of the interaction strength. They showed how to 
solve the infinite-dimensional DMFT model, illustrated 
how to determine the order-disorder transition temper- 
ature for a checkerboard (and incommensurate) charge- 
density-wave phase, showed how to find the free energy 
(including a first study of phase separation), examined 
properties of the spin-one-half model, and calculated the 
/-particle spectral function. 



as charge-density-wave 


order 


van Doneenl |l991a 


1992: van Donsen and Vollhardt 


. 119901 IFreericks 


199333; 


Gruber et all 


2001) and phase separa- 


tion llFreericks et all fl999 


iFreericks and Lemanski. 



Kimba ll problem of th e 
tion IjChung and Freericka 
was the problem of 



The original Falicov- 
met al- insulator transi- 
Il998ft was solved, as 
cla ssical intermediate va- 
lence (|Chung and Freericka l2000h [both using the 
spin-one-half gene r alizat ion (jBrandt et all Il990t 
IFreericks and Zlatid Il998ftl. The "Mott-like" metal- 
insulator transition l)van Dongen and Leinun d. Il997|) 
and the non-fermi-liauid behavior~ i|sT"er'aZ.l Il992ft were 
also investigated. Dynamical properties and trans- 
port have been determined ran ging fr om th e charge 
susceptibili ty ijFreericks and Millerl l2000l IShvaika , 
2000L 120011) . o the optical con ductivity llMoeller etal 



1992), to the Raman re spo nse llDevereaux et all 2003; 



2001a(), to an evaluation 


of the /-spectral func- 


tion (Brandt and Urbanek 


119921 ISi et all 11992 


Zlatic et all 12001). Finally, the static suscepti- 
bility for spontaneous polarization was also deter- 


mined iPortengen et all 


1996albl: ISi et al 


. 119921 


Subrahmanvam and Barma. 


19881 IZlatic et all 


2001). 



These solutions have allowed the FK model to be 
applied to a number of different experimental sys- 
tems i^n^in^_^ojn_^alence^hBiige-transition materi- 
als ijZlatic and Freericka l2001albl) like YbInCu 4 and 
EuNi2(Sii_ 2; Ge 2 ;)2 to materials that can be doped 
through a metal-insulator transi tion like Ta, T N [used as 
a barrier in Josephson junctions iFreericks et all [2001b, 
l2002dlMiiier and Freericka l200l[) ]. to Raman scattering 
in material s on the insulating side of the me tal-insulator 
transition l)Freericks and Devereauxl l2001b|) like FeSi or 
SmBg. The model, and some straightforward mod- 
ifications appropriate for double-ex change, have been 
used to describe the CMR materials llAUub and Alascic , 



| 1996lll997tlLetfulov and Freericka l200lt iRamakrishnaii , 
12003ft . 
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Generalizations of the FK model to the static Hol- 
stein model were first carried out by iMillis et al\ l)l995t 
1996) and also applied to the CMR materials. Later 
more fundamental properties relati ng to the transi- 
tion t emperature for the harmonic l)Blawid and Millisl 
2000; 



(of strength U) which represents the local Coulomb in- 
teraction when itinerant and localized electrons occupy 
the same lattice site (we could make U depend on a or 
•q, but this complicates the formulas and is not normally 
needed). The fourth term is the // Coulomb interaction 



Ciuchi and de PasauaTel 11999ft and anharmonic energy of strength UlL (which can be chosen to depend 



cases (Frcericks et al. 



2000 ) and relating to th e gap ra- 
tio for the har monic lIBlawid and Millis, 200l|) and an- 
harmonic cases (jFreericksandZlaticl 1200 laj) were com- 
pleted. Modifications to examine diluted magnetic semi- 
conductors have a l so ap peared l)Chattopadhvav et all 
l200lt iHwang et al\ . |2002|) . A new approach to DMFT, 
which allows the correlated hopping Falicov-Kimball 
model to be solved has also b een completed re- 
cently lISchilleTt 119991 IShvaikaL 12002ft . 



B. Hamiltonian and Symmetries 

The Falicov-Kimball model is the simplest 
model of correlated elect r ons. The original ver- 
sion ( Fa licov and Kimball Il969l) involved spin-one- 
half-electrons. Here, we will generalize to the case of 
an arbitrary degeneracy of the itinerant and localized 
electrons. The general Hamiltonian is then 



2s+l 



2S+1 



n 



-*EE c *W + E E E fnfU- 

(ij) ct=1 i rj=l 

2s+l 2S+1 

^EE E 4<wi/*i 

i a— 1 rj—l 
2S+1 

E/ E/ Uf t Jjr fjjjfirifir)' fir]' 
i rjrj' —1 

2s+l 

5MS^E E TO ^ C L C ^ 
i a— 1 
2S+1 

9/Mi3^E E m vfl V fiV 
i 77=1 



(1) 



The symbols c\ a (ci a ) denote the itinerant electron cre- 
ation (annihilation) operators at site i in state a (the 
index a takes 2s + 1 values). Similarly, the symbols fj 
ifiri) denote the localized electron creation (annihilation) 
operators at site i in state ij (the index r\ takes 25*+ 1 val- 
ues). Customarily, we identify the index a and r\ with the 
z-component of spin, but the index could denote other 
quantum numbers in more general cases. The first term is 
the kinetic energy (hopping) of the conduction electrons 
(with t denoting the nearest-neighbor hopping integral); 
the summation is over nearest neighbor sites i and j (we 
count each pair twice to guarantee hermiticity) . The sec- 
ond term is the localized electron site energy (which we 
allow to depend on the index rj to include crystal-field ef- 
fects [without spin-orbit coupling for simplicity]); in most 
applications the site-energy is taken to be ^-independent . 
The third term is the Falicov-Kimball interaction term 



on rj if desired; the term with r\ = rj is unnecessary and 
can be absorbed into Ef v ). Finally, the fifth and sixth 
terms represent the magnetic energy due to the inter- 
action with an external magnetic field H, with /j,b the 
Bohr magneton, g (<?/) the respective Lande g- factors, 
and m a (m ri ) the z-component of spin for the respective 
states. Chemical potentials fi and /i/ are employed to ad- 
just the itinerant and localized electron concentrations, 
respectively (in cases where the localized particle is fixed 
independently of the itinerant electron concentration, the 
localized particle chemical potential [if can be absorbed 
into the site-energy E$\ in cases where the localized par- 
ticles are electrons, they share a common chemical po- 
tential with the conduction electrons \i — fif). 

The spinless case corresponds to the case where s = 
5 = and there is no // interaction term because of the 
Pauli exclusion principle. The original Falicov-Kimball 
model corresponds to the case where s = S = 1/2, with 
spin-one-half electrons for both itinerant and localized 
cases (and the limit XJ*> — > 00). 

The Hamiltonian in Eq. possesses a number of dif- 
ferent symmetries. The partial particle-hole symmetry 
holds on a bipartite lattice in no magnetic field (H = 0), 
where the lattice sites can be organized into two sublat- 
tices A and B, and the hopping integral only connects 
different sublattices. In this case, one performs a partial 
particle-hole symmetry transf ormation on either the itin - 
erant or localized electrons l)Kennedv and LiebL Il986|) . 
The transformation includes a phase factor of (—1) for 
electrons on the B sublattice. When the partial particle- 
hole transformation is applied to the itinerant electrons, 



(-1) 



p(i) 



(2) 



with p(i) — for i £ A and p(i) = 1 for i € B and h 
denoting the hole operators, then the Hamiltonian maps 
onto itself (when expressed in terms of the hole operators 
for the itinerant electrons), up to a numerical shift, with 
U — > — U, Ef n — > Ef v + U, and — > — fx. When applied 
to the localized electrons, 



(3) 



the Hamiltonian maps onto itself (when expressed in 
terms of the hole operators for the localized electrons), 
up to a numerical shift with U — > —U,fi^ fx + U, 

Efn - ~E fv - E v > Ki> - U "r and W -> -W- 
These particle-hole symmetries are particularly useful 

when Ef v — 0, , U^L does not depend on 77, and we 
work in the canonical formalism with fixed values of p e 
and pf, the total itinerant and localized electron densi- 
ties. Then, one can show that the ground-state energies 
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of 7i are simply related 

E g . s .(p e ,p f ,U) = E g . s .(2s + 1 - p e ,pf,-U ) 
= S g . s .( Pe ,25+l-p/,-C/) 
= £ 9 . s .(2s + l-p e ,2S+l- P/ ,C/), 

(4) 

(up to constant shifts or shifts proportional to p e or pj) 
and one can restrict the phase space to p e < s 4- \ and 

Pf <s + i 

When one or more of the Coulomb interactions are in- 
finite, t hen there are additional symm etries to the Hamil- 
tonian ijFreericks et all Il999t l2002bl . In the case where 
all Lfff = oo, then we are restricted to the subspace 
Pf < 1. This system is formally identical to the case 
of spinless localized electrons, and we will develop a full 
solution of this limit using DMFT below. The extra sym- 
metry is precisely that of Eq. Q , but now with S = ( 
regardless of the number of r\ states). Similarly, when 
both = oo and U = oo, then Eq. (QJ holds for 

s = S = as well (regardless of the number of a and 
r\ states). These infinite-t/ symmetries are also related 
to particle-hole symmetry, but now restricted to the low- 
est Hubbard band in the system, since all upper Hubbard 
bands are pushed out to infinite energy. 

The Falicov-Kimball model also possesses a local sym- 
metry, related to the localized particles. One can easily 
show that [H, fj n fin] = 0, implying that the local occu- 
pancy of the /-electrons is conserved. Indeed, this leads 
to a local U(l) symmetry, as the phase of the localized 
electrons can be rotated at will, without any effect on 
the Hamiltonian. Because of this local gauge symmetry, 
Elitzur's theorem requires that there can be no quantum- 
mechanical mixing of the .f-particle number at finite 
temp erature l)Elitzuitll975llSubrahmanvam and Barmal 
1988), hence the system can never develop a spontaneous 
hybridization (except possibly at T = 0). 

There are a number of different ways to provide a phys- 
ical interpreta tion of the Falicov-K i mball model. In the 
original idea IjFalicov and Kimball [1969J), we think of 
having itinerant and localized electrons, that can change 
their statistical occupancy as a function of temperature 
(maintaining a constant total number of electrons) . This 
is the interpretation that leads to a metal-insulator tran- 
sition due to the change in the occupancy of the different 
electronic levels, rather than via a change in the charac- 
ter of the electronic states themselves (the Mott-Hubbard 
appro ach) . Another interpretation l)Kennedv and Liebl 
1986), is to consider the localized particles as ions, which 
have an attractive interaction with the electrons. Then 
one can examine how the Pauli principle forces the sys- 
tem to minimize its energy by crystallizing into a pe- 
riodic arrangement of ions and electrons (as seen in 
nearly all condensed-matter systems at low tempera- 
ture) . Finally, we can map ont o a binary alloy prob- 
lem llFreericks and FalicovL ^990) , where the presence of 
an "ion" denotes the A species, and the absence of an 
"ion" denotes a B species, with U becoming the differ- 
ence in site energies for an electron on an A or a B site. 



In these latter two interpretations, the localized particle 
number is always a constant, and a canonical formalism 
is most appropriate. In the first interpretation, a grand 
canonical ensemble is the best approach, with a common 
chemical potential (/i = p,f) for the itinerant and local- 
ized electrons. 

In addition to the traditional Falicov-Kimball model, 
in which conduction electrons interact with a discrete set 
of classical variables (the localized electron-number op- 
erators), there is another class of static models that can 
be solved using the same ki nd of techniques — the static 
anharmonic Holstein model ( HolstH 119591: iMillis et all 
119951) . This is a model of classical phonons interacting 
with conduction electrons and can be viewed as replacing 
the discrete spin variable of the FK model by a continu- 
ous classical field. The phonon is an Einstein mode, with 
infinite mass (and hence zero frequency), but nonzero 
spring constant. One can add any form of local anhar- 
monic potential for the phonons into the system as well. 
The Hamiltonian becomes (in the spin one-half case for 
the conduction electrons, with one phonon mode per site) 

Ti-Hoi = -t c\ a c j<y + g ep ^2 Xi{c\ a c ic! - p ea ) 

(ij)o icr 

+ ^ K ^2 X i + A» l X! X i + aan X! X i> ( 5 ) 

i i i 

where, for concreteness, we assumed a quartic phonon po- 
tential. The phonon coordinate at site i is Xi, g ep is the 
electron-phonon interaction strength (the so-called defor- 
mation potential), and the coefficients j3 an and a an mea- 
sure the strength of the (anharmonic) cubic and quartic 
contributions to the local phonon potential. Note that 
the phonon couples to the fluctuations in the local elec- 
tronic charge (rather than the total charge) . This makes 
no difference for a harmonic system, where the shift in 
the phonon coordinate can always be absorbed, but it 
does make a difference for the anharmonic case, where 
such shifts cannot be transformed away. The particle- 
hole symmetry of this model is similar to that of the 
discrete Falicov-Kimball model, described above, except 
the particle-hole transformation on the phonon coordi- 
nate requires us to send Xi — > — X4. Hence, the presence 
of a cubic contribution to the phonon potential nn ^ , 
break s the particle-hole symmetry of the system l)Hirschl 
1993); in this case the phase diagram is not symmet- 
ric about half filling for the electrons. We will discuss 
some results of the static Holstein model, but we will not 
discuss any further extensions (such as including dou- 
ble exchange for colossal magnetoresistance materials or 
including interactions with classical spins to describe di- 
luted magnetic semiconductors). 



C. Outline of the Review 

The Falicov-Kimball model (and the static Holstein 
model) become exactly solvable in the limit of infinite 
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spatial dimensions (or equivalently when the coordina- 
tion number of the lattice becomes large). This occurs 
because both the self energy and the (relevant) irre- 
ducible two-particle vertices are local. The procedure in- 
volves a mapping of the infinite-dimensional lattice prob- 
lem onto a single-site impurity problem in the presence of 
a time-dependent (dynamical) mean field. The path in- 
tegral for the partition function can be evaluated exactly 
via the so-called static approximation (in an arbitrary 
time-dependent field). Hence the problem is reduced 
to one of "quadratures" to determine the correct self- 
consistent dynamical mean field for the quantum system. 
One can next employ the Baym-Kadanoff conserving ap- 
proximation approach to determine the self energies and 
the irreducible charge vertices (both static and dynamic). 
Armed with these quantities, one can calculate essentially 
all many-body correlation functions imaginable, ranging 
from static charge-density order to a dynamical Raman 
response. Finally, one can also calculate the properties of 
the /-electron spectral function, and with that, one can 
calculate the susceptibility for spontaneous hybridization 
formation. The value of the Falicov-Kimball model lies 
in the fact that all of these many-body properties can be 
determined exactly, and thereby form a useful benchmark 
for the properties of correlated electronic systems. 

In Section II, we review the formalism that develops 
the exact solution for all of these different properties em- 
ploying DMFT. Our attempt is to provide all details of 
the most important derivations, and summarizing formu- 
lae for some of the more complicated results, which are 
treated fully in the literature. We believe that this review 
provides a useful starting point for interested researchers 
to then understand that literature. Section III presents 
a summary of the results for a number of different prop- 
erties of the model, concentrating mainly on the spinless 
and spin-one- half cases. In Section IV, we provide a num- 
ber of examples where the Falicov-Kimball model can be 
applied to model real materials, concentrating mainly on 
valence-change systems like YbInCu4. We discuss a num- 
ber of interesting new directions in Section V, followed 
by our conclusions in Section VI. 



II. FORMALISM 

A. Limit of Infinite Spatial Dimensions 

In 1989, Metzner and Vollhardt demonstrated that 
the many-b ody problem simplified i n the limit of large 
dimensions ijMetzner and Vollhardtl IT989(k equivalently, 
this observation could be noted to be a simplifica- 
tion when the coordination number Z on a lattice be- 
comes large. Such ideas find their origin in the jus- 
tification of the inverse coordination number 1/Z as 
the small parameter governing th e convergence of the 
coher e nt potential approximat i on jSchwartz and S iggia. 
11972ft . | Metzner and Voll hardlJ l)l989|) introduced an im- 



portant scaling of the hopping matrix element 

t = t*/2Vd = t*/V2Z (6) 

where d is the spatial dimension. In the limit where 
d — > oo, the hopping to nearest neighbors vanishes, but 
the coordination number becomes infinite — this is the 
only scaling that produces a nontrivial electronic den- 
sity of states in the large dimensional limit. Since there 
is a noninteracting "band", one can observe the effects 
of the competition of kinetic-energy derealization, with 
potential energy localization, which forms the crux of 
the many-body problem. Hence, this limit provides an 
example where the many-body problem can be solved ex- 
actly, and these solutions can be analyzed for correlated- 
electron behavior. 

Indeed, the central limit theorem shows that the nonin- 
teracting dBnsitv_of_s_tate3_^^ lattice phy P (e) 



tcractmg density oi states on a nybercut 
satisfies ilMetzner and Vollhardtl, 1 1989|) 



Ph VP {e) 



■exp(-eVO. 



(7) 



which follows from the fact that the band structure is 
a sum of cosines, which are distributed between —1 and 
1 for a "general" wave-vector in the Brillouin zone (f2„ c 
is the volume of the unit cell, which we normally take 
to be equal to one). Adding together d cosines will pro- 
duce a sum that typically grows like \fd, which is why 
the hopping is chosen to scale like \j\fd. The central 
limit theorem then states that the distribution of these 
energies is in a Gaussian [an alternate derivation rely- 
ing on tight-binding Green's function s and the proper- 
ties of B essel functions can be found in lMuller-Hartmannl 
l|l989af) ]. Another common lattice that is examined is the 
infinite-coordination Bethe lattice (which can be thought 
of as the interior of a larg e Cayley tree). T he noninter- 
acting density of states is ijEconomoul Il983ft 



PBethe(c) 



1 



2nt* 2 n 



(8) 



where we used the number of neighbors Z = Ad and the 
scaling in Eq. ©. 

The foundation for DMFT comes from two facts: first 
the self energy is a local quantity (possesses temporal 
but not spatial fluctuations) and second it is a functional 
of the local interacting Green's function. These obser- 
vations hold for any "impurity" model as well, where 
the self energy can be extracted by a functional deriva- 
tive of the Luttinger-Ward skelet on expansion for the 
self- e nergy generating functional ijLuttinger and Wardl 
1960). Hence, a solution of the impurity problem pro- 
vides the functional relationship between the Green's 
function and the self energy. A second relationship is 
found from Dyson's equation, which expresses the lo- 
cal Green's function as a summation of the momentum- 
dependent Green's functions over all momenta in the 
Brillouin zone. Since the self energy has no momen- 
tum dependence, this relation is a simple integral relation 
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(called the Hilbert transformation) of the noninteracting 
density of states. Combining these two ideas together, 
in a self-consistent fashion provides the basic strategy of 
DMFT. 
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FIG. 1 Skeleton expansion for the itinerant electron self en- 
ergy Eij through fourth order. The wide solid lines denote 
itinerant electron Green's functions and the thin solid lines 
denote localized electron Green's functions; the dotted lines 
denote the Falicov-Kimball interaction U. The series is iden- 
tical to that of the Hubbard model, except we must restrict 
the localized electron propagator to be diagonal in real space 
(which reduces the number of off-diagonal diagrams signifi- 
cantly). The first diagrams in the second and third rows are 
the only diagrams that contribute when i ^ j for the Falicov- 
Kimball model. 



For the Falicov-Kimball model, we need to establish 
these two facts. The locality of the self energy is es- 
tablished most directly from an examination of the per- 
turbation series, where one can show nonlocal self ener- 
gies are smaller by powers of l/Vd. The skeleton ex- 
pansion for the self energy (determined by the functional 
derivative of the Luttinger-Ward self-energy generating 
functional with respect to G), appears in Fig. ^ through 
fourth order in U. This expansion is identical to the ex- 
pansion for the Hubbard model, except we explicitly note 
the localized and itinerant Green's functions graphically. 
Since the localized electron propagator is local (i.e. no 
off-diagonal spatial components), many of the diagrams 
in Fig. n are also purely local. The only nonlocal di- 
agrams through fourth order are the first diagrams in 
the second and third rows. If we suppose i and j cor- 



respond to nearest neighbors, then we immediately con- 
clude that the diagram has three factors of Gy oc 1/ y/d? 
(each l/Vd factor comes from Uj). Summing over all 
of the 2d nearest neighbors, still produces a result that 
scales like l/Vd in the large dimensional limit, which van- 
ishes. A simi lar argument can be extended to all nonlo - 
cal diagrams ijBrandt and Mielscblll989tlMetznerl ll991h 
Hence, the self energy is local in the infinite-dimensional 
limit. The functional dependence on the local Green's 
function then follows from the skeleton expansion, re- 
stricted to the local self energy (i = j). 



B. Single-Particle Properties (Itinerant Electrons) 

We start our analysis by deriving the set of equations 
satisfied by the single-particle lattice Green's function 
for the itinerant electrons defined by the time-ordered 
product 

(9) 

where j3 = 1/T is the inverse temperature, < r < (3 
is the imaginary time, Zl is the lattice partition func- 
tion, ./V is the total itinerant electron number, Nf is the 
total localized electron number, and %- denotes imagi- 
nary time ordering (earlier times to the right). The time 
dependence of the electrons is 



C«t(t) 



_ p r(H-^N) 



c icr (0)e 



-t{H-hN) 



(10) 



and the trace is over all of the itinerant and localized 
electronic states. It is convenient to introduce a path- 
integral formulation using Grassman variables V>i<r(r) and 
V'io- ( r ) f° r the itinerant electrons at site-i with spin a. 
Using the Grassman form for the coherent states, then 
produces the following path integral 



G ija {r) = -±Tr f e-^f-^ 
J P^V^(T)Vv(0)e- s 



(11) 



for the Green's function, where TLj denotes the /- 
electron-only piece of the Hamiltonian [corresponding to 
the second, fourth and sixth terms in Eq. QJ]. The lat- 
tice action associated with the Hamiltonian in Eq. i jXj l 
is 



2s+l 



ij a=l 

c d 



t* 



x W(t'), (12) 



with Nfi the total number of /-electrons at site i. Since 
the Grassman variables are antiperiodic on the interval 
< t < /3, we can expand them in Fourier modes, 
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indexed by the fermionic Matsubara frequencies iu) n = 
inT(2n+l): 

OO 

i/>Ur(r)=T e- iu " T ^ ia {iw n ) (13) 



71— — OO 



and 



tpia{r)=T ^2 e lWnr ipi a (iu>„). 



(14) 



71— — OO 



In terms of this new set of Grassman variables, the lattice 
action becomes S L —T Y]°°_ ~, Si/, with 



2s+l 

EE 

ij a=l 
t 



2Vd 



{—iu>n - fi - g\iBHm a + UNfi)6ij 
lpi tr (iWn)'<Pj<7( iu3 n) (15) 



and the Fourier component of the local Green's function 
becomes 



dre^ T G iia {T) 



Zl 



x j P)/iD^ M (jw n )i/) !a (jw rl )e T E n ' s n '. 

(16) 

Now we are ready to begin the derivation of the DMFT 
equations for the Green's functions. We start with the 



many-body- variant of the cavity method ijGeorges et all 
1996; iGruber et oil l200lTl . where we separate the path 
integral into pieces that involve site-i only and all other 
terms. The action, is then broken into three pieces: (i) 
the local piece at site i, S n (i, i); (ii) the piece that couples 
to site i, S n (i,j); and the piece that does not involve site 
i at all, 5 n (cavity) (called the cavity piece). Hence, 



S n = S n (i,i) + S n (i,j) + S„(cavity) (17) 



with 



2s+l 



S n (i,i) = (-ioj n — M ~ g^BHm a + UNfj) 

X IpiaiilV^lpiaii^n), (18) 



(19) 



and 



5„ (cavity) 



S n (i,i) - S n (i,j). 



(20) 



Let Z cav ity and (— ) C avity denote the partition function 
and path-integral average associated with the cavity ac- 
tion S'cavity = Tj2 n 5n (cavity) (the path- integral average 
for the cavity is divided by the cavity partition function) . 
Then the Green's function can be written as 



J 



2s+l 



X / exp|T ^ ^2 E -^jj$i<T'(iu n ,)'lpj cr >(iu) n ,) + $jo-> (iUn'tyia' fan')]} > ■ (21) 

\ n' = -oo jXij^O = ! I cavity 

A simple power-counti ng argument, shows th at only the second moment of the exponential factor in the cavity average 
contributes as d — > oo ijGeorges et all Il996|) , which then becomes 

oo 2s+l 

x exp[T ^2 E {iu n > + M + g^BHm^ - UNfj- (ioy j}^/ (toy )Vw (joy)]. (22) 

n'=— oo a' — 1 

with \i a (iuj n ) the function that results from the second-moment average that is called the dynamical mean field. On a 
Bethe lattice, one finds Xi a (ioj n ) = t* 2 Ga a (iuj n ) , while on a general lattice, one finds \i a (iuj n ) — J2jk tijtik[Gjka(i^n) — 
Gji„(iu) n )Giko{iu n ) i 'Gii a (iuj n )} ijGeorees et all . \l996^i . But the exact form is not important for deriving the DMFT 
equations. Instead we simply need to note that after integrating over the cavity, we find an impurity path integral 
for the Green's function (after defining the impurity partition function via Zi mp = Zl/ 'Cavity) 

s-i , ■ s d\n Zi mp 

G iia (tLJ n ) - 



d\ ia {iuj n y 
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„ oo 2s+l 

Z imp = Tr f e- p Wt-*t"f) V$iDil>iexp[T £ ^ {iaj„ + /i + gp B Hrn a — U Nfi - X ia (iuj n )}-ipi a (iuj n )il; ia (iuj n )]. 

(23) 



n— — oo (7— 1 



r 



The impurity partition function is easy to calculate, be- 
cause the effective action is quadratic in the Grassman 
variables. Defining Z 0cr (/Lt) by 



iuj n + \i - A,, 



(24) 



n— — oo 



[where we used the notation \ na — A^ (ito„) and adjusted 
the prefactor to give the noninteracting result] allows us 
to write the partition function as 



-'imp 



Trfexpl-ftHfi-UfNfi)] 

2s+l 

Y[ Zoain + gp B Hm a - UN fl ). (25) 



[Zqv is the generating functional for the U = impurity 
problem and satisfies the relation Zg a = DetGj^T, with 
GqJ defined below in frequency space in Eq. It 
is cumbersome to write out the trace over the fermionic 
states in Eq. (|25ll f or the gene r al cas e. The spin-one-half 
case appears in lBrandt et al Here we consider 

only the strong- interaction limit, where — > oo, so 
that there is no double occupancy of the /-electrons. In 
this case, we find 

2s+l 



Zimp = J], Z a{^ + g^BHm^) 



<T=1 

2S+1 

7,=1 

2s+l 

I [ Z 0a {p + gp B Hm a - U). 



(26) 



Evaluating the derivative yields 

w 



Gna 



with 



Wo 



iu> n + fj, + gn B Hm a - X na 
Wl 

iuj n + fi + g\x B Hm a - \ na - U ' 



2s+l 

= Yl Zoaip + gp, B Hm a )/Z 

imp i 



(27) 



(28) 



and wi = 1 — wq. The weight w\ equals the average /- 
electron concentration pf. If we relaxed the restriction 
— * oo, then we would have additional terms cor- 
responding to Wi, with 1 < i < 25 + 1. Defining the 
effective medium (or bare Green's function) via 

[Goaiiun)}^ 1 = iu n + A 1 + g^BHm a - A nCT , (29) 



allows us to re-express Eq. <|2"?|) as 

wq wi 



G nr7 



(30) 



which is a form that often appears in the literature. 
Dyson's equation for the impurity self energy is 

S„ CT = T, a (iuj n ) = [G 0CT (iw„)] _1 - [G a (iu n )]~ x . (31) 

This impurity self energy is equated with the local self 
energy of the lattice (since they satisfy the same skele- 
ton expansion with respect to the local Green's function). 
We can calculate the local Green's function directly from 
this self energy by performing a spatial Fourier trans- 
form of the momentum-dependent Green's function on 
the lattice: 



Gna — ^ G W cr(k) 



E 



i 



iw n + pt, + gp B Hm a - S„ CT - e k 



dep(e) 



1 



iuj n + n + gp B Hm a - E n(T - e 



(32) 



where ek is the noninteracting bandstructurc. This re- 
lation is called the Hilbert transformation of the nonin- 
teracting density of states. Equating the Green's func- 
tion in Eq. J3"2"fl to that in Eq. l|3Uj) is the self con- 
sistency relation of DMFT. Substituting Eq. (|3*TJl into 
Eq. H30(l to eliminate the bare Green's function, pro- 
du ces a quadratic equation for the self energy, solved 
bv lBrandt and Mielschl (|l989l fl990> 



U- 



Gn 



± 



U- 



Gn 



4w 



U 



G r , 



(33) 

which is the exact summation of the skeleton expansion 
for the self energy in terms of the interacting Green's 
function (the sign in front of the square root is chosen 
to preserve analyticity of the self energy). Note that 
because w± is a complicated functional of G n and £„, 
Eq. (|33|l actually corresponds to a highly nonlinear func- 
tional relation between the Green's function and the self 
energy. 

There are two independent strategies that one can 
employ to calculate the FK model Green's functions. 
The original method of Brandt and Mielsch is to 
substitute Eq. (|3^1) into Eq. which produces 

a tra nscendental equation for Gna i n the complex 
plane (jBrandt and Mielschl Il989l Il990|) (for fixed w x ). 
This equation can be solved by a one-dimensional com- 
plex root-finding technique (typically Newton's method 
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or Mueller's method is employed), as long as one pays at- 
tention to maintaining the analyticity of the self energy 
by choosing the proper sign for the square root [for large 
U, the sign c hanges at a critic al value of the Matsub- 
ara frequency l)Freerickslll99 3a)]. The ot her technique i s 
the iterative technique first introduced bv l.Iarrelil jl992), 
which is the most commonly used method. One starts the 
algorithm either with S ntT = or with it chosen appro- 
priately from an earlier calculation. Evaluating Eq. I|32l) 
for G na then allows one to calculate Go from Eq. I|31l) . If 
we work at fixed chemical potential and fixed Ef, then we 
must determine wq and w\ from Eq. 128fl : this step is not 
necessary if wq and wi are fixed in a canonical ensemble. 
Eq. (|30|l is employed to find the new Green's function, 
and Eq. Ij31|l is used to extract the new self energy. The 
algorithm then iterates to convergence by starting with 
this new self energy. Typically, the algorithm converges 
to eight decimal points in less than 100 iterations, but for 
some regions of parameter space the equations can cither 
converge very slowly, or not converge at all. Convergence 
can be accelerated by averaging the last iteration with 
the new result for determining the new self energy, but 
there are regions of parameter space where the iterative 
technique does not appear to converge. 

Once wq, wi, and p are known from the imaginary-axis 
calculation, one can employ the analytic continuation of 
Eqs. (with iu n — > to + i0 + ) to calculate G a (u) 

on the real axis. In general, the convergence is slower on 
the real axis, than on the imaginary axis, with the spec- 
tral weight slowest to converge near correlation-induced 
"band edges" . A stringent consistency test of this tech- 
nique is a comparison of the Green's function calculated 
directly on the imaginary axis to that found from the 
spectral formula 



G a (z) = duj- 
J > 



to - 



M+' 1 



with 



and 



A a (uj) = / dep(e)A a (e,w), 



(34) 



(35) 



A <T (e,u) = Im — 

7T LO - 



1 



p + gfiBHm a - So-(w) — e • 



i0+ 
(36) 

[the infinitesimal + is needed only when ImS CT (a;) = 0]. 
Note that this definition of the interacting DOS has the 
chemical potential located at u> = 0. 

In zero magnetic field, one need perform these calcu- 
lations for one a state only, since all G n(T 's are equal, 
but in a magnetic field, one must perform 2s + 1 parallel 
calculations to determine the G nCT 's. 

Our derivation of the single-particle Green's functions 
has followed a path-integral approach throughout. One 
could have used an equation of moti on approach instead . 
This technique has been reviewed in lZlatic" et all (|2001[) . 



T he formalism for th e static Holstein model is sim- 
ilar ijMillis et all . 0.996). Taking 77 to be a continuous 
variable (x) and H = 0, we find 



.00 2s+l 

Zi m p = I dx Y\_ Z 0rT (p - g ep x) 
J — 00 1 



x exp(-(3[-Kx 2 + f] an x 3 + a an x i }), (37) 



and 



2s+l 



dx 



w(x) 



XLOjl ~\~ p QepX A? l( j 



(38) 



W ^ = IT 2 CT (M - 9epX) 
a=l 

x exp(-/3[-Kx 2 + (3 an x 3 + a an x 4 ])/Z lm p-(-i9) 

These equations l|37H39|l and Eqs. JsHJ and <EH are all 

that are needed to determine the Green's functions by 
using the iterative algorithm. 

The final single-particle quantity of interest is the 
Helmholtz free energy per lattice site Fffelm.- We re- 
strict our discussion to the case of zero magnetic field 
H = and to Ef v = Ej (no 77 dependence to Ef). 
There are two equivalent ways to calculate the free en- 
ergy. The original method uses the noninteracting func- 
tional form for the free energy, with the interacting den- 
sity of state s A(e) replacing the noninteracting d ensity of 
states p(e) l)PlischkeL H97I iRamirez et a/.l . ll970|) 

F He im. (lattice) - (2s + 1) J def{e)A(e)(e + p) + E f n f 

+ (2s + l)T J de{/(e)ln/(e) 

+ [l-/(e)]ln[l-/(e)]}^(e) 



+ T[w\ In wi + wq In wq — wi ln(25 + 1)] , 



(40) 



where /(e) = 1/[1 + exp(/3e)] is the Fermi-Dirac distribu- 
tion function. This form has the itinerant electron energy 
(plus interactions) and the localized electron energy on 
the first line [the shift by p is needed because A(e) is 
defined to have e = lie at the chemical potential p], 
the itinerant electron entropy on the second and third 
lines, and the localized electron entropy on the fourth 
line. Note that the total energy does not need the stan- 
dard many-body correction to remove double counting of 
the interaction becaus e the localized particles co mmute 
with the Hamiltonian l|Fetter and WaleckaL fl971^ . 

Th e Brandt-Mielsch app roach is differ- 
ent ijBrandt and Mielschl Il99lj) and is based on 
the equality of the impurity and the lattice Luttinger- 
Ward self-energy ge nerating fun ctionals <&. A general 
conserving analysis ijBavml 1 1962") shows that the lattice 
free energy satisfies 

FHeim. (lattice) = T^i att - E n(T G no - 
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tW dep(e) In 

71 (7 

fip e + flfWl 



1 



p + gp B Hm a - E ncr - e 
(41) 



whil e the impurity (or atom ic) free energy satis- 
fies l)Brandt and MielschLll99l 

F H eim. (impurity) = T$ imp - T^Yl na G na 

TUT 

+ Hp e + HfW\. (42) 

It a 

The first two terms on the right hand side of Eqs. I|41|) 
and 142(1 are equal, so we immediately learn that 

F H eim. (lattice) = -T\nZ imp J dep(e) 



x ln[(iuj n + + gp B Hm a 



e)G na ] + pp e + pfWi 
(43) 



Many of these transitions are continuous (second-order) 
transitions, that can be described by the divergence of 
a static susceptibility at the transition temperature T c . 
Thus, it is useful to examine how one can calculate differ- 
ent susceptibilities within the FK model. In this section, 
we will examine the charge susceptibility for arbitrary 
spin, and then will consider the spin and superconducting 
susceptibilities for the spin-one-half model. In addition, 
we will examine how one can perform an ordered phase 
calculation when such an ord ered phase exists. Our dis- 
cussi o n follows closely that oflBrandt and Mielschl l|l989t 
1990); iFreericks and Zlatid l|l998|) and we consider only 
the case of vanishing external magnetic field H = and 
the case when Ef v has no rj dependence. 

We begin with the static itinerant-electron charge sus- 
ceptibility in real space defined by (our extra factor of 
2s + 1 makes the normalization simpler) 



since F He im. (impurity) = -Tin Z imp + pp e + p f wi. The 
equival ence of Eas. fiDl and 11131 has been explicitly 
shown ijShvaika and Freer icksi 1200 25. When calculating 
these terms numerically, one needs to use caution to en- 
sure that sufficient Matsubara frequencies are employed 
to guarantee convergence of the summation in Eq. (1431) . 



C. Static Charge, Spin, or Superconducting Order 

The FK model undergoes a number of different phase 
transitions as a function of the parameters of the system. 



X CC (R, - R, 



1 



2s + 1 
Tr c/ 



f> (e- 
dr[Tr c/ — 



<(tK(0)) 



Tr 



(e-P n nf) 



(44) 



where Zl is the lattice partition function, 



Ea S Jl 1( 4 C «T and n i( T ) 



exp[r(W 



pN)]n1 (0) exp[— t(H — pN)]. If we imagine intro- 
ducing a symmetry breaking field — hin c t to the 
Hamiltonian, then we can evaluate the susceptibility as 
a derivative with respect to this field (in the limit where 
h = Q) 



T s^S^ dG jja {iuj n ) _ T f-v. \ dG kia( iuj n) 

n a kl 

Since = [iu n + p + h k - E kka (iui n )]5ki - t% t /2Vd, we find 

X (R» - Rj) = 2s + i 1^1^ [~ Gij<7{iUn)Gjia{iu) n ) + 2^2^!^ G jk(J (iuj n )G kja (iuj n ) dG t iuo \ ^ 

n (7 k a' m kk<r m ) i 

, (Ml 



where we used the chain rule to relate the derivative of where we have defined 



the local self energy with respect to the field to a deriva- 
tive with respect to G times a derivative of G with respect 
to the field (recall, the self energy is a functional of the 
local Green's function). It is easy to verify that both 
Gij a Gji a and dGjj a /dhi are independent of a (indeed, 
this is why we set the external magnetic field to zero) . If 
we now perform a spatial Fourier transform of Eq. (|46ll , 
we find Dyson's equation 

X°M = XT (q) ~TJ2 Xn° (q)IT m X^(q), (47) 



,ccO 



V 



dhj 



(q) = ~y X] 2s + 1 ^Gijoi^GjiviiWr, 



x giq-(Ri-Rj) 

1 1 _ ^ dE CT (iw„) 



pCC 

nm ~ T2s + 1 ^ dG a ,Uu} m Y 

crcr' 



(48) 



V is the number of lattice sites and x cc (q) 
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T^ n Xn C (q)- The fact that we have taken a Fourier 
transform, implies that we are considering a periodic lat- 
tice here (like the hypercubic lattice); we will discuss 
below where these results are applicable to the Bcthe 
lattice. Note that it seems like we have made an as- 
sumption that the irreducible vertex is local. Indeed, 
the vertex for the lattice is not local, because the sec- 
ond functional derivative of the lattice Luttinger-Ward 
functional with respect to G does have nonlocal contri- 
butions as d — ► oo. These nonlocal corrections are only 
for a set of measure zero of q-values llGeorges et aZJ . fl996t 
iHettler et all l2000t IZlatic and Horvatid Il990|) . and one 
can safely replace the vertex by its local piece within any 
momentum summations, which is why Eqs. Ij47(l and Ij48(l 
are correct. 

The local irreducible vertex function T can be deter- 
mined by taking the relevant derivatives of the skele- 
ton expansion for the impurity self energy in Eq. (|33fl . 
The self energy depends explicitly on G n and implicitly 
through W\ . It is because w\ has G dependence that the 
vertex function differs from that of the coherent potential 
approximation (where the derivative of w\ with respect 
to G would be zero) . The irreducible vertex becomes 



-pec 

nm - ( 2s 



(7(7 K. 



( dYi na 



\ dwi ) G ^ \dGraa' , 



\8G n J Tnn 

dw\ \ 1 



Substituting the irreducible vertex into the Dyson equa- 
tion [Eq. 1071)], then yields 



cc ( \ = ccO ( x 1 ~ ( 9S " CTl /9wi)G„ gl 7(q) 

X»W Xn W 1 + x ccO {q){d X n<7i/dGn<7i)wi 

with the function 7(q) defined by 

dwi 



, (50) 



7(q)-E^«(q)E 



dG n 



(51) 



and we have chosen a particular spin state <j\ to eval- 
uate the derivatives in Eq. (|5*0jl since they do not de- 
pend on a\. Multiplying Eq. (|5*0f by ^2 a (dwi/dG n cj) 
and summing over n yields an equation for 7(q). Defin- 
ing Z na = iuj n + fi - \ na = G~l + Y, na and noting that 
dwx/dGn = ^2 m [dwi/dZ m ][dZ m /dG n ], allows us to re- 
place dwi/dG n by 



(52) 





l ~ G n (fir) 






dG n r2 


i dwi 1 9£ m 
1 l^m dZ m \ dwi j 





(49) and solve the equation for 7(q) to yield 



7(qj 



E Ttff dw 1 /dZ na [l - Gl(d^ n /dG n ) Wl \/[l + G K 77»(q) - G? n (dY. n /dG n ) Wl ] 
1 - E n . dw l /dZ n<T G n t ln {ii){dT ln /dw 1 ) G J[l + G„77„(q) - Gl(dX n /dG n ) % 



(53) 



with rj n (q) defined by 
?7 n (q) = G n 



Gl X5i c0 (q) 



(54) 



and we have dropped the explicit a dependence for G 
and S. The full charge-density- wave susceptibility then 
follows 

cc(a] = _ T V [l-7(q)(9S ra /^ 1 ) G JG^ 
X W ^l + G n r) n (q)-Gl{dE n /dG n ) Wl ' 

(55) 

The derivatives needed in Eqs. (|53|) and (|55|) can be 

determined straightforwardly: 



2s+l 

E 



dwi 



{2s+l)w 1 {l-w 1 )UG 2 n 



dZna (1 + G„S„)[1 + G„(£„ - U)] ' 



(56) 



1-G 



2 ( 9Yi n 



\dG ri 



(l + G n £„)[l + G n (E n -E/)] 



l + G n (2H n -U) 



(57) 



and 



G 



2 fd^r. 



UGl 



l + G„(2£„-£7) 



(58) 



The final expression for the susceptibility and for 7(q) 
appears in Table [I] 

The mixed c/ susceptibilities can be calculated by 
ta king derivatives of Wi with resp ect to hi as shown 
in iBrandt and Mielscbl l)l989T) and iFreericks and Zlatid 
(1998). We won't repeat the details here, just the end re- 
sult in Table [I] The calculation of the // susceptibilities 
is similar. By recognizing that we could have calculated 
the cf susceptibility by adding a local /-electron chem- 
ical potential and taking the derivative of the itinerant 
electron concentration with respect to the local field, we 
can derive a relation between the cf and // susceptibil- 
ities. These results are also summarized in the table. 

In addition to charge susceptibilities, we also can cal- 
culate spin and pair-field susceptibilities for s > 0. We 
consider the spin-one-half case in detail. The spin sus- 
ceptibility vertex simplifies, since the off-diagonal terms 
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now cancel, and one finds a relatively simple result. 
The mixed cf spin susceptibility vanishes (because the 
Green's function depends only on the total /-electron 
concentration). The // spin susceptibility is difficult to 
determine in general, but it assumes a Curie form for 
q = 0. The pair-field susceptibility is determined by em- 
ploying a Nambu-Gor'kov formalism and taking the limit 
where the pair-field vanishes. There is an off-diagonal 
dynamical mean field analogous to A„ and w\ depends 
quadratically on this off-diagonal field. Hence, in the 
normal state, the irreducible pair-field vertex is diagonal 
in the Matsubara frequencies, just like the spin vertex. 
This means the susceptibility is easy to calculate as a 
function of the bare pair-field susceptibility 

xS B °(q) = ~ E G« J -t(H)%(-iwn)e* 1 ' (R '- B * ) . 

(59) 

The result appears in Table [I] 

The cc charge susceptibility diverges whenever 7(q) di- 
verges, which happens when the denominator in Tabled 
vanishes. Since this denominator is identical for the cc, 
cf, and // susceptibilities, all three diverge at the same 
transition temperature as w e would expect. This yields 
the s ame result as found in iBrandt and Mie lsch (1989, 
1990) for the spinless case, except there is an additional 
factor of 2s + 1 multiplying the sum in the denominator 
arising from the 2s + 1 derivatives of w\, which are all 
equal. This factor "essentially" increases T c by 2s + 1 
over that found in the spinless case. The existence of 
a T c is easy to establish, since the summation in the de- 
nominator goes to zero like 1/T 4 for large T and diverges 
like C/T for small T. If C > 0, then there is a transition. 

The q-dependence of the charge susceptibility comes 
entirely from rj n (q) and hence from the bare susceptibil- 
ity 



} (q) = -^G n (k+q)G„(k) 



1 



de- 



Pie) 



x F, 



y/l - X 2 J-oo iu> n + ii - £„ - e 
iuj n + /J, — E„ — Xe 



(60) 



VI ~X 2 

where all of the q-dependence can be summarized 



in a single param e ter X( a) = lim<i-»c>o J2i=i cos Qi/d 
j|Miiller-Hartmannl Il989a|) and we use F ao (z) — 
J dep(e)/(z — e) to denote the Hilbert transform. The 
results for Y° c0 (q) and y n (q) simpl ify for three gen- 
eral cases (IBrandt and MielscbL Il989|) : X = — 1 , which 
corresponds to the "checkerboard" zone-boundary point 
Q = (tt, 7r, it, ...); X = 1, which corresponds to the uni- 
form zone-center point q = 0; and X — 0, which corre- 
sponds to a general momentum vector in the Brillouin 
zone [since the value of the cosine will look like a random 
number for a general wave vector and the summation will 
grow like Vd, implying X — ► 0; only a set of measure zero 
of momenta have X(q) ^ 0]. The results for Xn C °(q) an d 
T] n (q) appear in Table ITT1 for the hypercubic lattice. Both 
the uniform (X = 1) and the "checkerboard" (X = — 1) 
susceptibilities can be defined for the Bethe lattice, but 
there does not seem to be any simple way to extend the 
definition to all X. Indeed, higher-period ordered phases 
on the Bethe lattice see m to have first - order (discontin- 
uous) phase transitions l|Gruber et all l200lj) so such a 
generalization is not needed. 

We will find that near half filling, the X = — 1 charge 
susceptibility diverges and far away from half filling the 
X = 1 charge susceptibility diverges. The X = sus- 
ceptibility never diverges at finite T. This is also true 
for the spin and pair-field susceptibilities. They are al- 
ways finite at finite T. One can easily understand why 
the spin susceptibility does not diverge — it arises simply 
from the fact that the spins are independent of each other 
and do not interact. Similarly, one can understand why 
the pair- field suscep t ibility does not diverge — Anderson's 
theorem ijAndersonl Il959t iBergmann and Raineit Il974j) 
states that one cannot have superconductivity with a 
static electron-electron interaction; the interaction must 
be dynamic. 

Since the (X = —1) "checkerboard" charge sus- 
ceptib ility diverges, we can also examine the ordered 
state l)Brandt and Mielschl . fl990). In this case, we have 
a charge-density wave, with different electronic densities 
on each of the two sublattices of the bipartite lattice. 
Hence both G^ CT (G^) and (S^) are different on the 
two sublattices. Evaluating the momentum-dependent 
Green's functions yields 



G A ' B (q) = iuj n + ^ + gfi B Hm a - + e q 

which can be summed over q to yield 

qA = f dc _ P( € ) iuj n + \i + gfi B Hm a - Y% a ^ = f de _ P( € ) iu} ™ + M + gp B Hm a - S^ CT 



(iuj„ + [i + gp B Hm a - E^ CT )(iw„ +fi + gfi B Hm a - E^ CT ). (62) 
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TABLE I Static charge, spin, and pair-field susceptibilities for the Falicov-Kimball model. The charge susceptibility is given 
for the general case, the spin susceptibility for spin-one-half (and only q=0 for the // spin susceptibility), and the pair-field 
susceptibility only for spin-one-half in the cc channel. Note that the cf charge susceptibility is equal to 7(q). 



Charge 



[l+G„(2S„-t7)- 7 (q)l/]G 2 

|l+G„(2E„-(7)]G„r l „(q) + (l+G„E„)[l + G„(E„-C7)] 



X cc (q) = -TEr= 
X c/ (q) =7(q) 

_ E^ = -^(2s+l)n, 1 (l-u, 1 )UGl/{ll+G n (2-E n -U)]G n y n (ci) + (l+G n j: n )[l+G n (-S n -U)]} 

~~ l-E~ = _ 00 (2«+l)»l(l-™l)^G^„( q )/(l+G„E„)[l+G„(E„-!7)]{[l+G„(2E„-i7)]G„ J 7, 1 (q) + (l+G„E„)[l+G„(E Il -!7)]} 

X W ~~ l-E~ = _ =0 (2s+l)»l(l-»l)^Ga t7 „(q)/(l+G„E„)[l+G„(E„-!7)]{[l+G„(2E„-i7)]G„7 J ,,(q) + (l+G„E,,)[l + G„(E„-i7)]} 



Spin 



^/nl-.TV 00 G^[l+G„(2S„-[/)] 

X \H) 1 Z^n=-oo [l+G„(2E„-(7)]G„r)„(q) + (l+G„E„)[l+G„(E„-(7)] 

x' c/ (q) = o 

X ,// (q = 0) =vn/2T 



Pair-field r(q)=TE: 



Cn°(q){: 



u 1 (l-m 1 )i7 2 |l + G, 1 S, 1 | 2 |l + G„(S, 1 -i7)r 



|l+G„(E„-[l- mi ]!7)| 2 [(l-™ 1 )|l+G„(E„-C/)| 2 + u , 1 |l+G„E„| 2 ] 



TABLE II Values of Xn C °(q) and ?7„(q) for the special X 
points 1, 0, and —1 on the hypercubic lattice. 



*(q) 


Xn C °(q) 




»7n(q) 


—i 


~G n /(iuj n + fj,- 


E„) 














i 


2[1 — (iuj n + (J,— 


E n )Gr n ] 


~ gT 217 



The algorithm to solve for the Green's functions is mod- 
ified to the following: (i) start with a guess for 
and ( or se t both to zero); (ii) evaluate Eq. ljtj2l) to 
find G^a and then determine G^ tT (i(jj n ) from Eq. H31|) 
evaluated on the A sublattice [S^ CT = {G^ CT (icj„)} _1 — 
(G^ CT ) -1 ]; (iii) determine Wq and w^ 1 from the appropri- 
ate A-sublattice generalization of Eq. (|28[1 ; (iv) evaluate 
Eq. (301 on the A-sublattice to find G^ a and Eq. (j2U to 
find S^ CT ; (v) now find G^ CT from Eq. JH2> and G^ a (iuj n ) 
from Eq. i|31[) : (vi) determine W(f and from Eq. I|28|) 
on the B-sublattice; and (vii) evaluate Eq. (|3U[) on the 
B-sublattice to find Gf la and Eq. to find S^ ff . Now 
repeat (ii-vii) until convergence is reached. We will need 
to adjust Ef — fif until (wf + wf )/2 — w\. The cal- 
culations are then finished and the order parameter is 
(wf — Wi)/2. A similar genera lization can be employed 
for t he static Holstein model l|Ciuchi and de Pasaualel 
1999). Since the skeleton expansion for the self energy 
in terms of the local Green's function is unknown for the 
static Holstein model (and hence it is not obvious how 
to determine T) , this is the most direct way to search for 
the ordered checkerboard phase. 



D. Dynamical Charge Susceptibility 

We now examine the dynamical charge suscep tibil- 
ity llFreericks and Milled. 1200(1 IShvaikai 1200(1 l200lh . de- 
fined by 



x cc (q,m) 



1 / dre iviT - V e q(R ^ R ^ 
23 + 1 Jo V nh 3 

(e-^<(rK(0)} 



Tr c/ 
Tr c/ - 



-Tr c/ - 



(63) 



with iv\ = 2iirlT the bosonic Matsubara frequency. Once 
again we assume we are in zero magnetic field H = and 
the /-electron site-energy is ^-independent Ef n = Ef. 
Our analysis follows the path-integral approach-one can 
also employ a strong -coupling perturbatio n theory to de- 
rive these formulae IjShvaikal 1200 0. 2001). The dynam- 
ical susceptibilities of the FK model have some subtle 
properties to them that arise from the fact that the lo- 
cal /-electron number is conserved. In particular, the 
isothermal susceptibilities (calculated by adding an ex- 
ternal field to the system, and determining how it mod- 
ifies the system in the limit where the field vanishes) 
and the so-called isolated (Kubo) susceptibilities [which 
assume the system starts in equilibrium at zero field, 
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then is removed from the thermal bath (isolated) and 
the field is turned on slowly] differ from each other (the 
response of the isolated system to the field is the iso- 
lated suscep tibility) due to the con s erved nature of the 
/-electrons i|ShvaikaL 1200 ll IWilcoxt Il968|) . In particu- 
lar, the isolated susceptibility vanishes for the // and 
mixed cf susceptibilities [due to the fact that the local 
/-electron number is conserved (i.e., [H,n{] = 0)], but 
is nonzero for the cc susceptibility. The divergence of 
the static charge susceptibility arises entirely from the 
coupling between the itinerant and localized electronic 
systems; the pure itinerant electron response (isolated 
susceptibility) never diverges. Hence the isothermal sus- 
ceptibility is discontinuous at zero frequency (i.e., it is 
not analytic). The continuous (analytic) susceptibility 
is the isolated susceptibility, and we will spend our time 
discussing it for the conduction electrons. 

If we express the susceptibility as a matrix in Matsub- 
ara frequency space, we find the following Dyson equa- 
tion: 

X cc (q.,iu m ,iu; n ;in) 
= x cca {^i^ m ]ivi)8 mn - T^x cc0 (q,^m; in) 

n' 

x T(iuj m ,iui n i;ivi)x cc {q.,iun' ,iu n ;ivi), (64) 

and the susceptibility is found by summing 
over the Matsubara frequencies x cc (it ^ v i) = 
T I2mnX cc (. c li'i' UJ m ,iLu n ;iv l ). Once again, the bare 
susceptibility depends only on the momentum parame- 
ter X, and takes the form 



which reduces to our static results for / = 0. The bare 
dynamical susceptibility simplifies in three cases that are 
summarized in Table ITTT1 Note that one needs to evalu- 



ate x with l'Hopital's rule whenever the denominator 
vanishes and we dropped the spin subscripts in the Ta- 
ble. The irreducible charge vertex is calculated for the 
impurity, and it satisfies 



r(iuj m ,iu! n ;ii>i) 



2s 



— V 

+ 1 z - 



1 8Y, a (iu m , iiv m+ i) 
T 8G a (iuj n ,ioj n+ i) 



(66) 



where we now have both a self energy and a Green's 
function that depend on two Matsubara frequencies, be- 
cause these functions are not time-translation-invariant 
in imaginary time. This occurs because we need to add a 
time-dependent charge field — oItx{t) c t( T ) c v( T ) 
to the action in order to evaluate the dynamic charge sus- 
ceptibility and this time-dependent field removes time- 
translation invariance from the system (it does not de- 
pend on the time difference of the arguments of the 
fermionic variables) . 

TABLE III Values of x cc °(X, iu) m ; in) for the special X 
points 1, 0, and —1 on the hypercubic lattice. 

*(q) X cc °{X,im m ;ivi) 

— 1 —{G m + G m +i)/(iuj m + iio m+ i + 2/i — E m — E m+ ;) 

—G m G m +i 

1 — {G m — Gm+l) I (Wl + S m — T, m + l) 



x {X, iu m ;m) 



l 



2s 



krr 



1 



(k)G m+iCT (k+q) 
P(e) 



F, 



VT^x^J-c 

ibJm+l + A* ' 



IU) 

E 



-Xe 



Vi -x 2 



E m e 

(65) 



It is not easy to perform calculations for Green's func- 
tions th at depend on two time varia bles. The original 
work by IBrandt and Urbanekl l)l992j) illustrates how to 
proceed. We start with the definition of an auxiliary 
Green's function 



97 X (r,r') 



Tr c T T / e -P Ho exp 


E s fo d TX(f)cUf)c,(?)_ 


c.(r)4(r')) 




e 0n ex p 


Iodf X (f) 


|2s+l 



(67) 



where Ho = — A* c \ c <? an d the time-dependence is 
with respect to Ho [the auxiliary time-dependent field 
x( T ) — Si x{i v i) exp(— iz//r) should not be confused with 
any susceptibility] . It is easy to show that this auxiliary 
Green's function is antiperiodic with respect to either r 
variable being increased by [3. Hence, we can perform a 
double Fourier transform to yield 



g* ux (iw m ,iu n ) 



= T [ dr [ T dr'e^ r 5 ° ua; (r,r')e~^" T ',(68) 

and this is the same Matsubara-frcqucncy dependence 
as in Eq. I|66|l . Substituting in the fermionic Grassman 
variables from Eq. I|13|l and restricting ourselves to x 
fields that satisfy xi} v o) = produces the following path- 
integral form for the auxiliary Green's function: 



g™ x {iw m ,iu n ) = - 
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x exp 



(69) 



where Z aux = [l+e^ fJ -] 2s+1 is the auxiliary partition func- 
tion for x(^o) = 0. We calculate the Green's function by 
adding an infinitesimal field Txmna^ncr'^ma and noting 
that g° ux (iu} m ,iu! n ) = d\nZ aux /d\ mncr . We will restrict 
our discussion to the case where only one Fourier com- 
ponent xfal) is nonzero. The matrix then has a nonzero 
diagonal, and one nonzero off-diagonal, whose elements 
are all equal to xfai)- Now the partition function is the 
determinant of the matrix that appears in the action of 
Eq. I|69|) , which assumes the simple form of the product of 
all diagonal elements (multiplied by a constant to assure 
the correct limiting behavior) 



gaux 



2s+l 



(70) 



When we add the extra field Xa to the action, the eigen- 
values of the matrix change only for m = n, m + I = n, 
and the a value of the %<r field. The former case is the di- 
agonal matrix element, which is simply shifted by Xmma- 
The shift in the latter case can be worked out via pertur- 
bation theory to lowest order in X- We find the mth and 
m + lth eigenvalues change to iuj m + ^i — Xm+imX{i v i) I 
and iuj m+ i + fj, + Xm+lrnXfai) / f ivi, respectively. Taking 
the relevant derivatives to calculate the auxiliary Green's 
function, then yields 



g% vx (iu m ,iu n ) = 



S„ 



Sm+lnXiiVl) 



lW m + (1 
1 



IV I 



1 



(71) 



As before, we define Go by adding a \ a field to the 
action in Eq. (|69[1 T\ ncr t[) na -t[) na - (and change the par- 
tition fu nction normalization accordingly). It is easy 
to show ijBrandt and Urbanekl Il992t IZlatic et all 1200 1|) 
that due to the fact that one can write the partition 
function of a path integral with a quadratic action as 
Z = det g , with g the corresponding Green's function, 
one discovers 



Go" 1 = [g a 



XI. 



(72) 



The full Green's function involves an additional trace 
over /-states. The "itinerant-electron" piece of the im- 
purity Hamiltonian equals — fx ^ CT cJ.Co- when w = and 
equals (U — /i) cJ.c CT when rJ = 1. It is a straightfor- 
ward exercise to then show that 

G a (iuj m ,iu> n ) = w G 0<7 (iu> m ,iu> n ) + wi[G^ - Ul]^ n , 

with wo given by Eq. I|28|l for H = and w\ = 1 — wq. 
Note that all of the Green's functions in Eq. I|73[> are 



matrices, and the —1 superscript denotes the inverse of 
the corresponding matrix. 

Eq. H73JI can be rearranged by multiplying on the left 
or the right by matrices like G -1 , Gq 1 , and Gq 1 — U to 
produce the following two matrix equations (with matrix 
indices and spin indices suppressed) 

G,7 2 - (U + G-^Go 1 + (1 - w^UG- 1 = 0, (74) 



Go 2 - G^i 11 + G_1 ) + (1 " wiWG- 1 



0. 



(75) 



Adding these two equations together and collecting terms 
results in 



1 



G„ 1 --(U + G- 1 ) 



1 



(U+G^f + il-w^UG- 1 = 0. 

(76) 

Now we substitute the matrix self energy E for Go, with 
the self energy defined by Dyson's equation 



Y,(iuj m ,iuj n ) = [G ] mn - [G~ 
to yield a matrix quadratic equation 



(77) 



2"^ 



U) 



2(2 Wl - 1)C/G _1 + G~ 



(78) 

that relates the self energy to the Green's function. In 
the case where only one Fourier component / 7^ of % is 
nonzero, both the self energy and the Green's function 
are nonzero on the diagonal m — n and on the diag- 
onal shifted by I units m + I = n. Substituting into 
the quadratic equation, and solving for the shifted diag- 
onal component of the self-energy [to first order in x(i v l)] 
yields the amazingly simple result 



E(iu m ,iuj m+ i) = G(iw m ,iu m+ i) 



-Tn+l 



G m — G 



m-\-l 



(79) 



after some tedious algebra. In Eq. 179JI . the symbols E m 
and G rn denote the diagonal components of the self en- 
ergy and the Green's function, respectively, which are 
equal to the results we already calculated for the self en- 
ergy and Green's function when xfal) = 0- 

We can now calculate the irreducible dynamical charge 
vertex from Eq. I|66|) , which yields 



1 XL 



T G m — G 



(80) 



The dynamical charge vertex is then a relatively sim- 
ple object for I 7^ [the static vertex is much more 
complicated, as can be inferred from Eqs. I|48|l. 149J1 . 
and (|52l)1. The original de rivation, based on the atomic 
limit (Shvaika L 1200(1 12001ft . produced a much more com- 
plicated looking result for the vertex, but some tedious 
alge bra shows that the two forms are indeed identi- 
cal l|Freericks and Milled. 12000). Using this result for the 
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vertex, we immediately derive the final form for the dy- 
namical charge susceptibility on the imaginary axis 



X 



X (X, iuj m ;ivi) 



+ X cc0 {X 1 iuj m \wi) 



'G m -G m+l 

Note, that if we evaluate the uniform, dynamical charge 
susceptibility q = (X = 1), we find 



X 



%x = i,i H ^ ) = -Tj2 



G m — G m +i 
ivi 



0, (82) 



which is what we expect, because the total conduction 
electron charge commutes with the Hamiltonian, and 
hence has no r-dependence, implying only the static re- 
sponse can be nonzero. Note that x cc °{X = 1) ^ 0, 
the vertex corrections arc needed to produce a vanishing 
total susceptibility. 

In general, one is interested in the dynamical charge 
response on the real axis. To find the real (dynamical) 
response, we need to perform an analytical continuation 
of Eq. (|81|l to the real axis. Since the isothermal sus- 
ceptibility is discontinuous at v\ — 0, we can perform 
the analytic continuation only for the isolated suscepti- 
bility. The procedure is straightforward: we divide the 
complex plane into regions where the Green's functions, 
self energies, and susceptibilities are all analytic, and we 
express the Matsubara summations as contour integrals 
over the poles of the Fermi-Dirac distribution function. 
Then, under the assumption that there are no poles in 
the integrands other than those determined by the fermi 
factors, we deform the contours until they are parallel to 
the real axis. After replacing any fermi factors of the form 
f(Lu+ii>i) by /(w), we can make the analytic continuation 
ivi — > v + i0 + . The final expressions are straightforward, 
but cumbersome. We summarize them for X = — 1 and 
X = in Tabic llVl Only the X — — 1 result is meaning- 
ful for the Bethe lattice. Note that one could just as 
easily calculate the incommensurate dynamical charge 
response, but the equations become significantly more 
complicated. We will calculate a similar result when we 
examine inelastic X-ray scattering below. 

One can test the accuracy of the analytic continua- 
tion by employing the spectral formula for the dynami- 
cal charge susceptibility and calculating the charge sus- 
ceptibility at each of the bosonic Matsubara frequencies. 
Comparing the spectral form with the result found di- 
rectly from Eq. I|81|l is a stringent self-consistency test. 
In most calculations, these two forms agree to at least 
one part in a thousand. 



E. Static and Dynamical Transport 

Some of the most important many-body correlation 
functions to determine are the correlation functions re- 
lated to transp ort propert i es via their corre sponding 
Kubo formulas IIGreenwoo l lT^lKubollT95l . Trans- 
port can also be solved exactly in the FK model, and we 



illustrate here how to determine the optical conductiv- 
ity, the thermopower, the thermal conductivity, and the 
response to inelastic light scattering. 

The Kubo formula relates the response function to the 
corresponding current-current correlation functions. Be- 
fore we discuss how to calculate such correlation func- 
tions, it makes sense to describe the different transport 
currents that we consider. Each current (except the heat 
current) takes the generic form 



2s+l 



ia(q) = H 7a ( k + q / 2 ) c k+ qCT Cko 



(83) 



with 7a (k) the corresponding current vertex function. 
The heat-current operator takes the form 



2s+l 



kk' 



(84) 

with 7q (k) the vertex that arises from the kinetic energy 
and 7q (k, k ) the vertex from the potential energy. For 
conventional charge, and thermal transport, we are in- 
terested only in the q — > (uniform) limit of the current 
operators, but the finite-q dependence is important for 
inelastic (Raman) scattering of X-ray light. The current 
vertex functions are summarized in Tabic Ivi 

The Dyson equation for any current-current correlation 
function takes the form shown in Fig. [21 which is similar 
to that given by Eq. (|64[) . but the bare and interacting 
susceptibilities are now different due to the correspond- 
ing 7 a factors. Note that there are two coupled equations 
illustrated in Figs. El (a) and (b); these equations differ by 
the number of j a factors in them (of course, there is only 
one equation when 7 a = 1 as we saw for the dynamical 
charge susceptibility). Note further, that one could eval- 
uate mixed current-current correlation functions where 
there is a j a on the left vertex and a 7b on the right ver- 
tex, but we don't describe that case in detail here. The 
irreducible vertex function T is the dynamical charge ver- 
tex of Eq. Q8U[) which has the full symmetry of the lat- 
tice. If the vertex factor 7 a does not have a projection 
onto the full symmetry of the lattice, then there are no 
vert ex corrections fr om the local dynamical charge ver- 
tex llKhuranallT99(^ . This occurs whenever 



7a (k + ^)G na (k)G n+la (k + q) = 0. 



(85) 



If Eq. H85|) is satisfied, then the corresponding current- 
current correlation function is given by the bare bubble 
depicted by the first diagram on the right-hand-side of 
Fig. El (a). 

One case where there are no vertex corrections is the 
optical conductivity, which is constructed from a q = 
current-current correlation function. We consider the 
case where the only effect of the magnetic field is from the 
Zeeman splitting of the states (which is accurate in large 
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TABLE IV Dynamical charge susceptibility on the real axis for X = — 1 and X = (the dynamical susceptibility vanishes for 
* = !)■ 

X Dynamical susceptibility \ cc (X,v) 



-1 J- r°° A, 1 f(,,\ [G(^)+G(^+^)]/[2^+2 M +^+S(^)-S(^+^)] 

1 2jr J -no UUJ \ J \ W ) l-[G(w) + G(w + i/)J[S(w)-i;(w + y)]/[2w + 2/i+y + E(tj)-E(tj + y)J[G(a;)-G(w + i/)] 

£(.._!_ ,.\ [G*(^) + G*(a> + i;)]/[2a J + 2 M + i/+S*0)-£*(^ + i;)] 

' l-[G*(w)+G*(w+i/)][E»(w)-E*(w+ 1 /)]/[2w+2;*+^+S*(w)-S»(w+i')][G»(a>)-G*(w+^)] 

r/r/-, A _ #(,,, _l ,,11 [G-(^)+G(^+^)]/[2^+2 M + 1 .+S*(^)-S(^+^)] 

U ^ ./ ^ l-[G»(w)+G(o>+«)][E»(w)-E(w+v)]/[2w+2M+«+S*(w)-E(w+v)][G*(a>)-G(w 

PI J. r 00 J, J ff, A G(u>)G(u>+v) r , , \ G*0)G*(a>+i;) 

u 2tt J-oc I J v") 1-G(^)G(^+v)\T,(u)S(u+v)\/[G(lo)-G(lo+v)\ * \ ' J 1-G* (w)G* (w) -S* (w+k)J/[G» (w)-G* (w+!/)J 

r r/ \ /•/ i \l G* (u;)G(u; + i^) 

-[Jl^J ~ A&> + ^JJ l-G*(u;)G(u; + y)[S*(w)-S(w + I /)]/[G*(u;)-G(u; + I /)] 



TABLE V Current vertex functions for use in Eqs. and 
1841 1. The symbols e ! a and e„ are the polarization vectors for 
the incident and outgoing light, respectively, and a denotes 
the lattice spacing. An overall factor depending on properties 
of the incident and outgoing light is neglected for the inelastic 
light scattering vertex. 

Current type vertex function y a (k) 

Charge 7™(k) = iea\7e(k)/h 

Heat Y Q (k) = i'aVe(k)[e(k) - fj]/h 

70(k, k') = iaU[X7e(k) + Ve(k')] W(k - k')/2h 
W(k) = exp(-ik ■ R rf ) / rf Vi,/V 

Inelastic 7i (k) = £ a/3 4^^^* 

light scattering 



dimensions, because the phase factors induced in the hop- 
ping matrix occurs only in one dimension) . We illustrate 
here how to calcu late the optical conductivity on the 
hypercubic lat tice ijMoeller et all Il992t iPruschke el all 
Il993albl fl995L The starting point is the current-current 
correlation function on the imaginary lattice, which is 
equal to the bare bubble diagram in Fig. (a) 

(J n k 

where we have chosen the 1 direction for the velocity op- 
erator and the n superscript denotes the particle-number 
(charge) current operator (note our energy unit t* equals 
1 so it is suppressed). The average of sin 2 (ki) (times a 
function of ek) is equal to 1/2 times the integral over e 
of the function times the density of states. Hence, the 
current-current correlation function becomes 

2 2s + 1 00 „oo 

a—l n— — 00 ^ 

1 

X 

iu n + [i + gfiBHrria - S ncr - e 




FIG. 2 Coupled Dyson equations for current- current corre- 
lation functions described by the vertex function j a - Panel 
(a) depicts the Dyson equation for the interacting correla- 
tion function, while panel (b) is the supplemental equation 
needed to solve for the correlation function. The symbol V 
stands for the local dynamical irreducible charge vertex given 
in Eq. I8UI . In situations where Eq. 1851 is satisfied, there 
are no charge vertex corrections, and the correlation function 
is simply given by the first diagram on the right hand side of 
panel (a). 



iuJn+i + M + g^BHm a - E n+ ia- - e 

Note that the optical conductivity is a 1/d correction. 
The analytic continuation is performed in the same fash- 
ion as was done for the dynamic charge susceptibil- 
ity. The optical conductivity is defined to be ct(lo) = 
lmx nn {v)/v, so the final result is 

o(v) — (Tq / dep(e) / dujA a (e, oj)A a (e, oj + v) 
a —i ^—00 J— 00 
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/M - f{u> + u) 



with (To = e 2 TT 2 /hda d 2 (which is approximately 4.3 x 
10 3 fi _1 -cin _1 in d = 3 with a w 3 x 10" 8 cm). Note that 
one factor of /i is needed to construct the dimensionless 
energy unit when multiplied by v in the denominator, 
which is why only one factor of h appears in ao- 

Eq. I|88|l can be simplified further by performing the 
integral over e. This is done by using Eq. I)36|) to rewrite 
both spectral functions as an imaginary part of a Green's 
function, and noting that the imaginary part of z is (z — 
z*)/2i. The integral over e can then be performed after 
expanding each integrand by partial fractions to yield 



o(y) 



a 2s+1 f 

y_ 



dujKe< 



G„{u)-G%(lu + v) 



^i/ + E ff (w)-E*(a; + i/) 
G a {u)-G a {u) + v) \ /(w) - /(w + i/) 



(89) 



On another lattice, the optical conductivity takes the 
form of Eq. (|88|l but with an extra factor of v 2 (e). En- 
forcing the optical sum rule (which relates the integral 
of the optical conductivity to the average kinet ic energy) 
i|Chattopadhvav et all 1200(1 iMaldaeuel Il977|) produces 
a differential equation for v 2 (e) 



-[ 

de l 



\e)p(e)} + ep(e) = 0, 



(90) 



with the boundary condition that v 2 (e) — > at the band 
edges (where appropriate). On the hypercubic lattice 
one needs to sum the optical conductivity over all d 
axis directions to be able to be put in the above form; 
solving the differential equation gives v 2 (e) = 1/2 as 
expected. For the infinite coordination Bethe lattice, 
one finds v 2 (e) = (4 — e 2 )/3. This agrees with the 



conje ctured result l)Chung and Freericksl Il998t IVelickvl 
1969), but it relies on enforcing the sum rule for the 
Bethe l attice, which has not been est ablished indepen- 
dently ijvan Dongen and Blumerl . l2002j) . 

If we take the limit v — > on the hypercubic lattice, 
we find 



&dc 



dw 



dm 

dio 



2s+l 

E r <»> ( 91 ) 

(T=l 



with the spin-dependent relaxation time equal to 
ImG CT (w) 



Im£ CT (u;) 

2Re{> + n + g f i B Hm a - £ CT (w)]G CT (a;)|92) 



In addition, one is often interested in thermal transport 
quantities such as the thermopower S and the electronic 
contribution to the thermal conductivity n e . These three 
quantities are usually expressed in terms of three different 
transport coefficients in, L12 = L21 and L22 as follows: 



Ode 



e 2 L u 



(93) 



S = 



\e\TLu' 



and 



h 2 
T 



L22 — 



L12L 
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(94) 



(95) 



with kg the Boltzmann constant. The individual trans- 
port coefficients are determined by the zero-frequency 
limit of the analytic continuation of the relevant polar- 
ization operators Lij = lim^^o Re[iLjj(z/)/V] with 



and 



13 



d T e WlT Ti, 



dre w,r Tt 



of 



of 







dTe w ' T Ti 



of 



(T T e- 


- m jn(r) Jn (0)) 




z L 


(T T e- 








(T T e- 


- m 3Q(r)j Q (0)) 


Z L 



(96) 



(97) 



(98) 



where the subscripts n and Q denote the number 
(charge) and heat currents respectively (and we sup- 
pressed the Cartesian vector indices). All of these cor- 
relation functions are determined by their correspond- 
ing bare bubbles, because there are no vertex correc- 
tions for any of them. Not e that our sign conven t ion fo r 
the thermopower is that of lAshcroft and Merminl l)l976tl . 
where the coefficient multiplying L^/Ln is positive for 
negatively charged carriers. This leads to the situa- 
tion where electron-like transport has a positive ther- 
mopower, and hole-like transport has a negative ther- 
mopower at low temperature. (The sign of the ther- 
mopower varies depending on what convention is used 
for the definition of S in terms of the tra n sport coeffi- 
cients.) A theorem bv I Jonson and Mahanl 1l98Cl Il990l) 
says that even in a correlated system, there is a sim- 
ple relation between these different transport coefficients 
[that they reproduce the so - called Mott noninteracting 
form iChester and Thelhinet Hofilty 



Ll} - e 2 



dus ( — 



dm 

dw 



2s+l 

E 

cr=l 



i+j-2 



(99) 



What is remarkable, is that in the case of the FK 
model, one can explicitly calculate the relevant correla- 
tion f unctions and verify di r ectly the Jons on-Mahan the- 
orem ijFreericks and Zlatid l2001bl l2002aD . The deriva- 
tion is too long to reproduce here. Note that we have 
taken the more modern definitions of the coeffi- 
cients here, which has one less power of T than the nor - 
malization used by lFreericks and Zlatid (|2001bl l2002a|) : 
iJonson and Mahanl Hl*98rill990|) . 
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The final transport property we determine is the scat- 
tering of inelastic light. When optical photons are used, 
this corresponds to conventional electronic Raman scat- 
tering: inelastic light scattering off of the charge exci- 
tations of the many-body system. When higher-energy 
(X-rays) are employed, one scatters the photon off of the 
momentum and frequency-dependent charge excitations 
of the system. The Raman scattering limit results in the 
limit q — > 0. For simplicity, we consider only the spinless 
case, with H = 0. 

Inelastic light scattering depends on the polarizations 
of the incident and outgoing light. As such, it provides 
some additional symmetry resolution over and above the 
elastic scattering of an optical conductivity measure- 
ment. There are traditionally three main symmetries 
considered in Raman scattering experiments: (i) A\ g 
which has the full symmetry of the lattice; (ii) B\ g which 
has a d-wave symmetry and (iii) Bi g which is another 
d-wave symmetry. Each symmetry is chosen by different 
polarizations for the incident and scattered light. Here 
we concentrate on nonresonant Raman scattering, where 
the vertex functions are not functions of the photon en- 
ergies. If we sum over the d pairs of polarizations, where 
e 1 = e and each vector points along each of the different 
Cartesian axes, then we have the A\ g sector. If we choose 
e 1 = (1, 1, 1, ...) and e° = (1,-1, 1, — 1, ...), then we have 
the Bi R sector. And if we choose e 1 = (1,0,1,0,...) 



1 



and 



Ml. 1.0. I .... . then we have the Bi g sec- 



tor. If we have just nearest-neighbor hopping, then the 
i?2 g response vanishes because 7s 2g = 0- Following the 
form given in Table |Vl we find 7^ (q) = — e(q) and 

7Bi g (q) =t* E^icosq^-Tp'/V^. 

A straightforward calculation, using Eq. 1)85(1 . shows 
that the B\ g response has no vertex corrections on the 
zone diagonal q — (q,q,q,q, ■••)• Hence, the zone- 
diagonal Bi g response is the bare bubble. The Ai g re- 
sponse (and the Bi g off of the zone diagonal), on the 
other hand, does have vertex corrections, and is more 
complicated. The calculation of each response function 
is straightforward, but tedious. One needs to first solve 
the coupled equations depicted in Fig.|21on the imaginary 
axis and then perform the analytic continuation as we 
have done previously for other response functions. The 
end result is quite long and is summarized in Tables IVTI 
and EIH 

It is interesting to take the limit of conventional Raman 
scattering with optical photons, where q — > (X — ► 1). 
In this case, one finds xo = — / dep(e)/[u> + p — — 
e] [lu + v + fi — £ (lu + v) — e] and the Bi g response simplifies 
to 



XB ls {v) = 7~ duj / de P( e ) 

J— 00 J — CO 



1 



lu + v + fi — S(w + v) — e 

1 1 



LU + /! — X(<x>) 



LU + (1 — E*(w) 



lu + fi — £*(u;) — e 
1 

lu + v + [i — Y,*{lu + v) — e 
1 

lu + v + fi — Tj(lu + v) — e 



(100) 



Now using the definition of the spectral function 
in Eq. (|36l) . and taking the imaginary part of 
Eq. i|100D produces the so-ca l led Shastry -Shraiman re- 
lation ijShastrv and Shraimanl Il990t Il99lj) which relates 
the imaginary part of the nonresonant Bi g Raman re- 
sponse function to the optical conductivity 



(101) 



which was first proved in IFreericks and Devereauxl 
fcOOlhl) . 

Another interesting limit is the q = (ir,Tv, it, ■■■) {X = 
— 1) limit, where we find Xo = Xo — 0> Xo = Xo/2, and 
Xo = Xo/2. In this case, the renormalization due to the 
dynamical charge vertex exactly cancels in the numera- 
tor and denominator, and one finds the same result for 
the inelastic X-ray scattering in the A\ g and B\ g sec- 
tors. Hence, a polarized measurement of the inelastic 
X-ray scattering at the zone boundary determines the 
relative importance of nonlocal charge fluctuations to the 
strongly correlated system. This is quantified by compar- 
ing the scattering in the two different symmetry sectors; 
the difference in the results is a measure of the nonlocal 
correlations. 

We have produced results for inelastic X-ray scatter- 
ing only along the zone diagonal. One could choose other 
directions as well. In general, we find that the symme- 
tries from different channels then mix, and both A\ g and 
Big channels are renormalized by the irreducible dynamic 
charge vertex. The formulas in that case are complicated 
and will not be presented here. 

One could also study resonant (or mixed) Raman re- 
sponse functions. The formalism is similar to the non- 
resonant case, except it is better to work directly on 
the real axis. These calculations have not yet been 
completed by anyone. But a simple power-counting 
analysis of the different symmetry sect ors can be per- 
formed l)Freericks and DevereauxLl2001bj) . What is found 
is that the A\ g sector has contributions from the nonres- 
onant, the mixed diagrams, and the resonant diagrams. 
The Bi g sector is either nonresonant or resonant, but the 
mixed contributions disappear as 1/d, The £?2 g response 
is purely resonant. 



F. Single- Particle Properties (Localized Electrons) 

We now study /-electron properties. Since the local 
/-electron number is conserved, one might believe that 
the localized particle dynamics are trivial. Indeed, this 
is true for the particle number, but the fermionic elec- 
tron degrees of freedom are nontrivial, even if they are 
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TABLE VI Response functions for inelastic X-ray scattering by a photon of momentum q = (q, q, q, ...) along the zone diagonal 
described by the parameter -X(q) = cosq for the Ai g and Bi g sectors. The symbols with xo appear in Table IViTl and * denotes 
complex conjugation. 

Scattering response function 

r., s , xo(^;^,^)+ g»i^:g < ( ^ ) ) [xo(^;X, v )xo(^;X, v )-x{f(^;X, v )} 
-UK") - n u + v )\ sy^-^+v, 

XBi B (q, = 37 JZ« d " {/(^)xo(w; x, ^ - /> + ^)xS(^; x, u) - [/H - /(w + «/)]xo(w; a, «/)} 



TABLE VII Symbols xo, Xo> X0i Xo, Xo> and xo that appear in Table IvTl 



Symbol 



X o(^;A» = - JZz dep(e) w+fi _l (ui) _ e ^_ x2 F^ 
Xo(u;X,u) = - dep(e) a)+M _J» M _ e F a 



+- 



tjj + v + /i — S ( uj + 1/ ) — X e 
+ i/+/i-i:(u + i')-Xe 

V^-x 2 

+v+/j,-'E(uj+v)-X£ 

V^x* 



-1 + 



oj + 1/+ — S ( w + ) 



V^x^ 



F a 



J + f + fJ-S(q) + l/)-X < 



+ 



j + p-£*(w)-e 



-1 + 



Xo(u>; V, „) = i X0 ( W; X, V ) - 1±2 dep(e) | ^-^p ^ji^f « ( 



+ 



+ 



W+ M -£(w)-e (1_X 2 )3/ 2 



V^ 2 



lj + v + — S ( a; + 1/ ) — X e 

Vi--^ 2 



uj-\-i/-\-fj, — 12(uj-\-v)—Xe 



0^f2 



+ 



-1 + 



OJ + + — S + ) 



F a 



o+v+^~T,(oj+u)-Xt 



)']} 



V^X~ 2 



j -\- v -\- fi — S(w + [ / ) — 



- } 



restricted to be local. As a first step in our analysis, 
we can imagine writing the statistical factor for the f- 
occupancy in terms of an effective /-level l(Czvcholi 
1999) 

Wl = 1 + (1 ° 2) 

which we have written explicitly for the spinless case, 
with suitable generalizations for higher-spin cases. No- 
tice that in many cases, either w\ changes as a function 
of T, or /j,f changes as a function of T, or both. Hence, 
E% typically has temperature dependence. If we view 
this energy as the centroid of the /-spectral function, wc 
can immediately learn about how the /-spectral function 
may change as a function of T. Of course one can only 
learn so much from a single number. 



What is more interesting is to evaluate the lo- 
cal /-electron Green ' s fun ction following the wor k of 
iBrandt and Urbanekl l|l992t) and lZlatic eTaj] (|200lh [see 
also iJanisI (|l99J) which examined the corresponding X- 
ray edge problem] . The local /-electron Green's function 
on the imaginary axis is defined by 



pn rp (e-^S(A)/,(T)/t(Q)) 

F v (t) = -Tr c/ , (103) 



for t > in the interaction representation [Zi mp is given 
by Eq. (|26l) with H = 0]. The impurity Hamiltonian in 
zero magnetic field is H lmp = — ^n c + J2„{Efv ~ ^f) n { + 
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Un c nf and the evolution operator satisfies 
rP rP 

S(A)=T T exp - dr dr' V 4(r)A a (r - r')c CT (r') 
Jo Jo a 

(104) 

with Ao-(t) = exp(-iw„r)A no -. The field \ na is 

the discrete function that satisfies Eqs. I|27l - I32|l . Since 
f v (r) = exp(rH lmp )/r ; (0) exp(-rH imp ), we find 



d 



Mr) = 



-(^- M/ )-[/V4(r) Ctr (r) 



Mr), 
(105) 



which can be integrated to yield 
Mr) = e-( E f«-"f> 

x T T cxp ^-C/^^ r dr' C t(r')c CT (r')^ /„(0). 

(106) 

If we define a time-dependent field 

X t (t',t'') = -C/0(t-t')£(t'-t"), (107) 

then the /-electron Green's function becomes 



1 



. e -(fi/,-W)r TrcTT ( e -^0 5(A _ Xr )) ; 



F v {r) = 

(108) 

with Hq = —[in c , because the trace over / states is re- 
stricted to rJ = for t > propagation [the operator 
/(0)/ f (0) projects onto n* = and all of the remaining 
operators are composed of /t/ pairs which do not change 
the total /-electron number] . 

The strategy for evaluating the trace over itinerant- 
electron states is the same as in Section lil.DI The trace 
involves operator averages that are not time-translation 
invariant because the Xt field is not a function of r' — t" 
only. We start with the auxiliary Green's function in 
Eq. (|69fl and note that the double Fourier transform of 
the Xt field [analogous to Eq. Q68[l] becomes 



xAm) 



T 



dr'e^'xriT 1 ) 



-^fl 



Ut 

' 



1^0; 
I = 0, 



(109) 



where our choice for normalizing the Fourier transform 
has an extra factor of T for convenience here [in other 
words Xrii^i) has dimensions of energy] and the bosonic 
Matusbara frequency arises from the difference of two 
fermionic Matsubara frequencies. Since Xrii^o) 7^ 0, we 
nowhaveZ a "(-Xr) = [l + exp(/3/Lt- Ut)] 2s+1 . The par- 
tition function can also be expressed as the determinant 
of the matrix that appears in the exponent in Eq. (|69ll . 
In the general case, we can write the partition function 
as an infinite product 

2s+l 

lLU n + fj, + XO 



: (-Xr) 



2e /3(M+Xo)/2 ~Q 



lUJr, 



(110) 



which is equal to the product of the diagonal elements of 
the corresponding matrix 



M mn = {iuj m + n)8 mn + Xr, 



(111) 



up to an overall normalization factor. Hence, the deter- 
minant of the matrix in Eq. fill]) is simply a product 
of the diagonal elements! This is surprising, since one 
would naively expect that the determinant depended on 
the higher Fourier modes of x- We can show this result 
explicitly, by checking that 



= T 



dre {i ' 



---')-exp{/ dr'[x(r') - Xo}} 
Jo 

(112) 

is the m'-th eigenvector, with eigenvalue iu> m i + /i + xo- 
The verification, is straightforward, provided one first 
writes Xm-n in terms of the integral over t, and notes 
that x( T ) — Xo times the last exponential factor in 
Eq. i|112|) is equal to the negative r derivative of the ex- 
ponential factor. The eigenvalue can then be determined 
by evaluating the integral by parts; we need to subtract 
the xo factor from the exponent in order for the bound- 
ary terms to vanish when one performs the integration 
by parts. 

Knowing an explicit form for the eigenvectors then al- 
lows us to determine the auxiliary Green's function in the 
general case. As before, we imagine adding an infinites- 



imal field Xr, 



X^mm Snn to the matrix M„ 



The 



auxiliary Green's function is determined by taking the 
logarithmic derivative of the partition function with re- 
spect to x: jr( w »o.O = dlnZ aux (~Xr)/dx- This 
derivative is easy to calculate, because it says we need 
only determine the shift in each of the eigenvalues for 
Z aux [— Xt) to first order in x- In the presence of x, it is 
easy to show that 

= iu m , + M + Xo + XttP*^ + 0(X 2 ), (H3) 
which leads then to 



°(iu n ,iu! m ) = ^ 



i(m')* ± (m') 
Om <Pn 



lU) m > + /! + xo 
m' — — 00 

r P 



x e 



Jo 

Jo df[x(f)- X o]„$o df[ x " (f)-xl] 



,(t-t') 



m — — 00 



iu m > + M + Xo 



(114) 



Note that if we choose x( T ) = xH v l) ex P( — i v i T )i with 
x{i v l) real, then evaluating Eq. (|114|l to first order 
in x{i y l) reproduces Eq. <|T1|> as it must. For our 
case, we are interested in the real x field of Eq. (|1U7|I . 
Substituting this field into Eq. I|114|) an d perform- 
ing some tedious algebra, then yields the results for 
the diagonal and off-diagonal auxiliary Green's function 
shown in Table IVIIII These results generalize those 
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TABLE VIII Auxiliary Green's function g% ux (iu) n , iu) m ) = T£oA nm +T(£o — l)B„ m for n / m and g% ux {iu> 
T(£ - l)D nn for n = m. Here £ = 1/[1 + exp([7r - f3fi)}. 

A 1 e i(u n +„)T-0» ( e -(^m+rir +e ~^ )t< ,(. M „+ri T - e CT ) t e -U" m -u n )T 

rLrn i^m— ^Ti)(*^T7i+M — U) (iuJ m -\-ii)(iuj n -\-[i) (ii J J m +fJ.)(ii J J n -\-fJ.— U) — LJ n ) (iLJ Tl + a0 — ^ti ) (i^Tn +^) 

r (i Um +c)T [i(Mm _ t)n)e i/T_ (i[|)[|i + )1 _ c , )e (i u „+|.)r| 

' i(w m -w„)(iwm+M-t r )(jw n +/i-t/) 
D _1 e -i(->m-".i)T _ _ e '"^^( e -^ m ^_ e (n~C/)^)_ e l3M( e ~(^ t?1 +^) T _ e -[/T ) i^„+ M -i(g )m -a,„) e ' 3 f-'''- , "'+^> T 

nm i(w m — w„)(iw n +^-(7) i(utn-Wn)(>"»+/<) (iw (ia>»+M) »Om -wn)(iw m +/i)(iLij„+/j) 

e -.( t . m - M „)r[ iLJn+)1 _ tf+1 ( 1Jm _ |Jn)e (i am +y-i;)T] 

i(u m -u n )(iu m +n-U)(iu> rt +ii-U) ^_ 

~iu)„+p-t/)T-l + e-t i ""+' J - t7 ) T ' l- e -('""+' J -' :f ) T +e-^' J (e< i ""+>') T -e i:fT ) . (iu,„ + M ) (,3-t) - 1 -e~f M+ +a * )t 



n -(^„+ M ~a)r-l + e'"^"+^- c/ > T , l_ e ('-n+^-P)T +e ji|.( e -[.»„+ l .) r _ e -g T ) | ( iQ , n+At) ( T _ /3) _ 1 _ e )3M-(^„ + ^)T 
-^nn 77~ j rn2 r 7T- i ,,irj,., _^,,_m I 



of lBrandt and Urbaneld il992T) off of half filling and cor- 
rect some tvpos in lZlatic et al\ l)200l|) . 

Once i?^"^ i s known, then we follow the same proce- 
dure as before to calculate the relevant averages. The 
object of interest is Tr c 7^-(exp(— f3H. )S(\ — Xt)), which 
can be viewed as the partition function corresponding 
to a particle moving in both the A and \ T fields. Since 
the partition function is the determinant of the inverse 
of the Green's function, and since the Green's function 
ga-(iuj n ,iuj m ) associated with the A — Xt field is related 
to the auxiliary Green's function via 



(115) 



we immediately conclude that 

Tr c T r (exp(-/?H )S(A-x T )) 

2s+l 

<T=1 

Substituting into Eq. (|108J) yields our final result for the 
/-electron Green's function 



Fr,(r) 



_ e -(%,-f/)'- 2 «(_ Xr ) 



Ul=i Det[«5 nm - g a a ux ^n, iw m )A mff ].(117) 



-'imp 



One needs to recall that Z aux (~ Xt) changes with r be- 
cause the Xt held depends on r. One can easily ver- 
ify in the case where Ef v has no rj dependence, and in 
the limit where r — > + , then g% vx — > 8 nm /(iuj n + n) 
and the /-electron Green's function can be evaluated ex- 
plicitly, producing F v (0 + ) = — w as it must. Similarly, 
for t — > (3~ one finds g% ux — > 5 nm /(iui n + \i — U) and 
F rl ifi-) = -w 1 . 
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FIG. 3 Contour c used in integrating the action in the 
Keldysh formalism. 



One is also interested in the /-electron Green's func- 
tion on the real axis. Unfortunately, there is no obvious 
way to analytically continue Eq. (|II7|) from the imagi- 
nary to the real axis. Instead, one is forced to re-evaluate 
the /-electron Green's function directly on the real axis. 
Since there is no time-translation invariance to the ef- 
fective fields, it is most convenient to employ a Keldysh 
formalism, even though the /-electron is in equilibrium. 
In the Keldysh formalism, we need to perform the path 
integral in the time domain over the contour illustrated 
in Fig.|21to determine the contour-ordered Green's func- 
tion. The Keldysh (greater) Green's function for real 
time t becomes 



F>(t) = -Tv cf - 



-^S c (A)/,(t)/t(0)> 



(118) 



with f v (t) = exp(itH im p)fri(0) eyLp(-itH imp ) and 



5 C (A) = T c exp 



dt Jdt' 4(*)A(i, i')<v(F) 



One can directly determine X(lu) on the real axis. Using 
this function one then finds 

1 f°° 

\(t,t') = / dwlmX(u)exp[-iu)(t-t')] 

n J-oo 

x [f{u)-e c {t-t')\ (120) 

where 6 c (i — t') = if i' is in front of i on the contour c 
and 1 if it is behind. 

Using an equation of motion, similar to Eq. I|1U5|) , leads 
us to introduce a time-dependent field 



Xt (i, i') = -iU0 c (t - t)S c (i- ¥). 



(121) 



The derivation of the Green's function follows as before 
resulting in the analogue of Eq. I|117|) 



F>(t) 



_t e ^ E i-n-^Y2, aux [—xt) 



■'imp 



2s+l 



x Y\ Det[S(t - t') 

+ j dt"gT' x {Vt")K{t",i% 



(122) 
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TABLE IX gZ ux (t,t') for different orderings of J, t, and t' 
along the contour c. The symbol £o satisfies £o = 1/[1 + 
exp(iUt — /3fJ,)]. 



svah(J3u)Tm{F > (t)} sia(wt)] 



(128) 



and 



£o exp[»>(t - 1')] 

fo exp[i(7t + — U)i — ijxt'\ 

foexp[i(A*-CO(t-f)] 

(|o - l)exp[i>(t-F)] 

(|o - 1) exp[-ii7t + i/ii - i(n - U)f) 

(to -l)exp[»0i -U)(t-F)] 



t<i <t' 
i<t<F 
t<¥ <t 
t<F<i 
F < t < i 
i ' < i < t 



1 f°° 

A fo (uj) = / dt[-smh(f3uj)Re{F>{t)}cos(ujt) 



except now, the determinant is of a continuous matrix 
operator. 

The function g" ux (i,t') can be found directly from its 
operator definition 



+ (cosh^u) + l)lm{F > (t)} sin(ajt)]. (129) 

At half filling, we have that A^°(w) = 0, so that 
Eqs. !gl§ and (|r2T)j> reduce to 

2 f°° 

A i n (uj) = — / dtRe{F>(t)}cos(u;t). (130) 
77 Jo 

This then determines the localized electron spectral func- 
tion on the real axis. 
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Ti c (T c e~ l3Ho S c (-xt)ca(t)ci(t')) G- Spontaneous Hybridization 



(123) 

noting that we order the times along the contour c and 
that Z aux (- X t) = [1 + exp(/3/i- iUt)] 2s+1 . The result is 
shown in Table ITXl for the six differ ent possible orderings 
of t, t , and t' along the contour c. iBrandt and Urbaneld 
(1992) show how to calculate the discrctized determinant 
in an efficient manner. If we use a quadrature rule 

I dbl{t) =J2WiI(ti) (124) 

with weights Wi for the discrete set of times {t{\ on the 
contour c, then the continuous determinant can be ap- 
proximated by the discrete determinant 



Det[W^ 



k 



**)WfcM*fc, *,•)}] (125) 



for each a (1/Ai c is the approximation to the delta func- 
tion on contour c with At c the width of the interval that 
includes the delta function; for a (midpoint) rectangular 
quadrature rule, one takes At c — 1/Wi). 

Now the Keldysh Green's function satisfies the follow- 
ing spectral formula 



due-™ *[/(«)- lJ^H 



(126) 



because it involves the product of the spectral function 
AUu)) with the relevant distribution function [f(io) — 1]. 
The (greater) Green's function obviously satisfies 



F>(t)=F>*(-t). 



(127) 



We break the spectral density into its even and odd 
pieces: A^iuj) = [A^lo) + At(-w)]/2 and Af°(cj) = 
[Af(oj) - Af(-u)}/2. Then we use Eq. ljrT7|) to show 
that 

1 f°° 

A f n e (uj) = / di[(cosh/?u; + l)Re{F>(i)}cos(u;i) 

n Jo 



Sham and coworkers proposed that correlations can 
lead to a spontaneous ferroelectricity via the creation 
of a d ynamically correlated hyb ridization in the FK 
model l|Portengen et all Il996albl) . Such an effect is 
necessarily a subtle one, because it is well known that 
Elitzur's theorem (Elitzur, 1975) applies to the FK 
model IjSubrahmanvam and Barmal Il988f) and no such 
spontaneous hybridization can occur at any finite tem- 
perature. Nume rical and analytic a l calculatio ns in 
finite-dimensions (Farkasovskv, 1997, 1999, 2002]) and 
in infinite-dimensions IjCzvcholi Il999[) indicated that 
such spontaneous hybridization did not occur in the FK 
model, but they did not cover all of the available param- 
eter space. An alternate way to test these ideas is to 
directly calculate the su sceptibility toward spontaneous 
hybridization formation l|Zlatic et all l200lh . 

Because the /-electron dynamics are local, only the 
local susceptibility toward spontaneous hybridization is 
relevant. We will restrict ourselves to the spinless case 
for simplicity. The local spontaneous hybridization sus- 
ceptibility is defined by 



Xhyb 



dr— 

o Z 



1 



imp 



x Tr c/ (% exp(-/?W imp )S(A)/(r)dt (r)d(0)/ t (0)) . 

(131) 

If we introduce the time-dependent field A(r, r') that cou- 
ples the itinerant electron at r' to its Hermitian conjugate 
at time r [in equilibrium, we have A(r, r') = A(r — r')], 
then we can express the susceptibility in terms of a func- 
tional derivative 



Xhyb — 



P e -{Ef-Hf—U)rga 

dr 



: (-Xr) $Det[l -g aux \] 



-'imp 



SX(t, 0) 



(132) 
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Noting that DetA = exp[Tr(ln A)], allows us to immedi- 
ately compute 



<5Det[l -g aux X] 
5X(t,0) 



-Vet[l-g aux X][{l-g aux Xy l g aux ]o T , 

. (133) 

where matrix multiplication is understood. Substituting 
Eq. (|133l) into Eq. l|132l) then immediately yields 



/ dr [ dr'F(r)(l - g aux \)^,g^, (134) 
Jo Jo 



Xhyb 



where one must recall that g aux must be recomputed for 
each t value, since it depends explicitly on Xt- Reex- 
pressing Eq. (|134fl in terms of Fourier transformed quan- 
tities produces 



Xhyb 



f drF(r)Tj2(l 9™* X) m \g aux ^ n ,r\ 

(135) 

where the partial Fourier transform of the auxiliary 
Green's function satisfies 

aux f ■ \ m \ r aux { - - \ iu^r fioc\ 

[iuj n ,T) =T 2_^g {iu) n ,iuj m )e m . (136) 
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Note that this final result for the spontaneous hybridiza- 
tion susceptibility in Eq. i|135|) requires the matrix in- 
verse of [1 — g aux X] which is straightforward to compute 
numerically for a finite truncation of the matrix. 



III. ANALYSIS OF SOLUTIONS 

A. Charge- Density- Wave Order and Phase Separation 

The first problem examined in the FK model with 
DMFT was the problem of ordering into a two-sublattice 
charge densit y wave (CDW) at half filli ng on a hyper- 
cubic lattice ijBrandt and Mielschi |T989). Because the 
hypercubic lattice is bipartite, we expect a transi t ion at 
finite temperature for all U l|Kennedv and LiehL Il986h 
and indeed this is true. Since the original work, a num- 
ber of ot her studies of CDW order and phase separation 
follow e d llBrandt and Mielschi Il990t 1 199 it Ivan Dongeri , 
19921 Ivan Dongen and yollhardtt Il990i iFreericks , 
1 993alli iFreericks^^IT i99ii"lFreericks and Lemafiski , 
200Gt iGruber et all l200ll iLetfulovl Il999|) on both the 
hypercubic and the Bethe lattices. We concentrate on 
the spinless FK model here. 

Since the DMFT is in the thermodynamic limit, we 
can determine the transition temperature for a continu- 
ous transition by simply finding the temperature at which 
the susceptibility diverges. Using the results of Table [I] 
we find that the susceptibility diverges whenever the de- 
nominator of 7(q) vanishes. In Fig. 0] we plot the T c 
for checkerboard CDW order (X = —1) and the spin- 
odal temperature for phase separation {X = 1) near half 
filling (a) and near the band edge (b). Note that the tran- 
sition temperature has a maximum that is about 1 /40th 



of the effective bandwidth (as At*) and at half filling for 
large U it behaves like T c sw t* 2 /AU [as expected for 
the equivalent (Ising) spin model] and for small U it ap- 
pears to grow like T c w U 2 \ In U \ [which is different from 
the expected exponentially growing behavi or that occurs 
in most nested systems at weak couplin g] ijvan Ponged . 
Il992t Ivan Dongen and Vollhardtl Il990j) . As we move 
away from half filling, the checkerboard phase is sup- 
pressed (disappearing at small and large U) and the 
phase separation is enhanced (especially at large U). 
Near the band edge [panel (b)], there is no checker- 
board instability, but the phase separation becomes even 
stronger. We see that the spinodal temperature for phase 
separation does not depend too strongly on the electron 
filling, and surprisingly, the curves cross as a function of 
U. 
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FIG. 4 T c for CDW order in the spinless FK model (a) near 
half filling and (b) near the band edge. The solid lines denote 
checkerboard order and the dashed lines denote the spinodal 
decomposition temperature for phase separation (the true 
first-order transition temperature is always higher than the 
spinodal temperature). The labels denote the electron filling 
p e [which runs from top to bottom for the dashed lines in 
panel (a) as 0.3, 0.35, 0.4, and 0.45, the spinodal temperature 
for 0.5 vanishes; in panel (b) the lines are for 0.05, 0.1, 0.15, 
and 0.2 — note how the curves cross as a function of U]; the 
ion filling is fixed at pf = 0.5. 

We can also examine the transition temperature for 
fixed U as a function of electron concentration (with 
Pf = 0.5 again) as shown in Fig. [5j For small U [panel 
(a)], we see that the region of the checkerboard phase 
increases as U increases, as does the region of phase sep- 
aration. For large U [panel (b)], we see that the T c curve 
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for the checkerboard phase develo ps a cusp at half fill- 
ing ijFreericks and Lemahskil l2000|) . and that the region 
of stability shrinks as U increases further. Phase separa- 
tion dominates for large U, as expected. 
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FIG. 5 T c for CDW order in the spinless FK model (a) for 
small U and (b) for large U. The solid lines denote checker- 
board order and the dashed lines denote the spinodal decom- 
position temperature for phase separation (the true first-order 
transition temperature is always higher than the spinodal 
temperature). The labels denote the value of U ; the ion filling 
is fixed &t pj — 0.5. 

It turns out, that for small U, the system also has in- 
stabilities to CDW phases with incommensurate values 
of X (i.e., X changes continuously with the electron con- 
centration). These re sults are summar ized with a pro- 
jected phase diagram IjFreeri cks. 1993a) in Fig.H3 which 
plots the regions of stability for different CDW phases 
as determined by the ordering wavevector of the initial 
ordered phase as T is lowered to the first (continuous) 
instability at T c (as determined by the divergence of the 
relevant susceptibility). This is an approximation of the 
zero-temperature phase diagram — the phase boundaries 
may change as one reduces the temperature from T c to 
zero, but we are not able to study the ordered phase of in- 
commensurate states. They may also change if there are 
first-order phase transitions. Note how the system segre- 
gates at large U for all fillings except p e = 0.5 where it 
is degenerate with the checkerboard phase. 

In order to understand the incommensurate order bet- 
ter, we plot T C (X) for different electron fillings at U = 0.5 
in Fig. [7| Near half filling, these curves are peaked at 
X = — 1, but as the system is doped sufficiently far from 
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FIG. 6 Projected phase diagram for the spinless FK model 
with pj — 0.5. The shaded region near half filling denotes the 
region of stability of the checkerboard phase (X — —1), the 
shaded region near the band edge denotes the region of sta- 
bility for the segregated phase (as determined by the spinodal 
decomposition temperature), and the white region shows the 
incommensurate phases at small U (the solid lines are lines of 
constant X). 



half filling, the curves develop a peak at an intermedi- 
ate value of X which evolves towards X — as T c — > 
(p e = 0.35 and 0.3). Then, if doped further from half fill- 
ing, T c rises again and the ordering wavevector evolves 
smoothly toward X = 1 (p e = 0.2 and 0.15). Note that 
the region of stability for positive (incommensurate) X 
is much smaller than that of negative (incommensurate) 
X. In the large U case, the maximum always appears 
at X = -1 or I = 1 and the incommensurate order 
disappears (see Fig. El . 

It is interes ting to also examine t he dynamical charge 
susce ptibility IjFreericks and Milled. 120001 IShvaikaL l2000l 
l200l|) . We don't expect the dynamical susceptibility to 
show any signs of the CDW instability because it is an 
isolated susceptibility, and the divergence of the static 
susceptibility arises solely from the coupling to the static 
electrons which produced an additional contribution to 
the (isothermal) static susceptibility. Indeed this is the 
case, as can be seen for the checkerboard (X = —1) dy- 
namical charge susceptibility at U — 4 shown in Fig. [5] 
We plot the imaginary part of the susceptibility only. 
Note how at high temperature there is a low-energy peak 
and a charge-transfer peak (centered at v w U = 4), 
but as the temperature is lowered, the low-energy peak 
rapidly disappears and the susceptibility has little fur- 
ther temperature dependence in the homogeneous phase 
(even though we have passed through T c which occurs 
at 0.0547). We also find that as we move from the zone 
corner toward the zone center, the higher-energy (charge- 
transfer peak) loses spectral weight but does not have 
much dispersion or peak narrowing, while the low-energy 
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FIG. 7 CDW transition temperature as a function of the 
ordering wavevector described by X. The numbers label 
the electron filling with the filling changing in steps of 0.05 
(U — 0.5 and p/ = 0.5). The results for p e — 0.25 are too low 
to be seen on this figure. 



peak shows an appreciable dispersion toward lower en- 
ergy and narrows as we approach the zone center (not 
shown). At the zone center, the dynamical charge sus- 
ceptibility vanishes. 
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FIG. 8 Imaginary part of the dynamical charge susceptibility 
at X = — 1 and [7 = 4. Note the low-energy spectral weight 
which rapidly disappears as T is lowered, and note further 
that there is no signal of the CDW order (which sets in at 
T c = 0.0547) in the isolated susceptibility. 



The results on the Bethe lattice are similar to those 
on the hypercubic lattice except for the incommensurate 
order. We will concentrate on discussing the Bethe lat- 
tice case for the remainder of this section. We begin by 
examining the ordered checkerboard phase at half fill- 
ing. In Fig. El we plot the electronic density of states for 
U = 1 above T c , and below T c in steps down to T = 0. 
Inset in the figure is a plot of the "normalized" order 
parameter — wf as a function of temperature. The 
order parameter takes a BCS-like form, being flat for low 
temperatures and having a square-root dependence as 
Tc is approached [ this behavior changes as U — > 0, see 
Ivan Doneenl (JT992)] . The density of states evolves with 



temperature in the ordered phase because of this tem- 
perature dependence in the order parameter. Each sub- 
lattice has a sharp peak in the density of states (which 
becomes singular at T = 0) as expected. Note the in- 
teresting behavior of the subgap states, which split into 
two subbands then shrink and disappear as T is reduced. 
In the limit as U — > oo, the bandwidth of the subbands 
goes to zero, and the densi ty of states become s a delta 
function on each sublattice {van Dongenlll992f) . 
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FIG. 9 Density of states above and below T c — 0.06614 in 
the spinless FK model at half filling on a Bethe lattice with 
U = l. The DOS is plotted for T = 0.07, T = 0.065, T = 0.06, 
T = 0.045, T = 0.03, and T = 0. Note that in the ordered 
phase, there are two DOS plots, one for each sublattice (solid 
and dotted lines); these DOS are mirror images of each other. 
The inset plots the "normalized" order parameter as a func- 
tion of temperature, which has a BCS-li ke shape. The exac t 
form is known for U — ► oo and U Ijvan Dongenl 1992), 
and one expects, for smaller U, that the shape will change 
dramatically since the order parameter becomes very small in 
the range 0.5T C < T < T c . 

In addition to the checkerboard phase, we can also 
examine the phase-separated (segregated) phase, where 
the itinerant electrons avoid the localized ions. If we 
take the limit U — > oo, where the phase separation 
is the strongest, we find the plot projected onto the 
wi = (1 — Pe)/2 plane (which corresponds to relative 
half filling of the lower Hubbard band) shown in Fig. 1101 
The solid line denotes the first-order phase transition line 
calculated by performing a Maxw ell construction on the 
free energy {Freericks ef aZlll999^ . The dotted line is the 
spinodal decomposition temperature calculated from the 
divergence of the susceptibility. Note how the two tem- 
peratures track well with one another and how they meet 
at the maximal T c as they must because that is the tem- 
perature where the first-order phase transition becomes 
second order. 

Finally, we exa mine higher-period ordered phases on 
the Bethe lattice IIGruber et a/1 |200T|) . Since the Bethe 
lattice does not have the same periodicity properties of 
regular lattices, one must construct higher-period or- 
dered phases with care. This is done by working directly 
in real space and repeating quasi-one-dimensional pat- 
terns of charge density waves according to the different 
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FIG. 10 Phase diagram for the first-order phase separation 
transition on a Bethe lattice with wi = (1— p e )/2 and U — > oo. 
The solid line is the first-order transition temperature and the 
dotted line is the spinodal temperature. 



levels of the tree that forms the Bethe lattice. A period- 
two phase is the conventional checkerboard phase with 
alternating charge densities on each level. There is a re- 
gion of parameter space where the period-three phase has 
been shown to be stable. This is illustrated in Fig. II II (a), 
which is a restricted phase diagram at T = that com- 
pares the ground-state energy of the segregated phase, 
the period-two phase, the period-three phase and the ho- 
mogeneous phase. Note the wide region of stability for 
the period-three phase (this region can shrink as addi- 
tional phases are added to the phase diagram; in partic- 
ular, we conjecture that the homogeneous phase stability 
will disappear as higher-period phases are added to the 
diagram). To examine in more detail, we plot the free 
energy of the homogeneous phase, the period-two phase, 
and the period three phase for U = 3, pt = 2/3, and 
/i = 0.84861 (p e w 0.332). Note how at low tempera- 
ture the period-three phase is lowest in energy, but how 
there must be a higher-period phase that intervenes be- 
tween the period-two and period three phases for the 
free energy to be continuous. We expect the phase tran- 
sition from the period-three phase to be first order. Note 
further how the period-two phase free energy smoothly 
joins the homogeneous-phase free energy, as expected for 
a continuous (second-order) phase transition. 



B. Mott-like Metal-Insulator Transitions 



a »0.5 

c 

o 



The FK model is not a fermi liquid <|Si et all Il992j) 

whenever pf is not equal to or 1 (when pf = or 
1, the FK model is a noninteracting fermi gas). This 
is because the /-electrons appear like disorder scatterers 
(with an annealed averaging rather than the conventional 
quenched averaging). These scatterers always produce a 
finite lifetime at the fermi surface, so rigorously quasipar- 
ticles do not exist and the system is a non-fermi-liquid. 
If U is small, one can still view the system as a "dirty 
fermi liquid" , but the system rapidly changes character 
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FIG. 11 (a) Restricted phase diagram on the Bethe lattice 
for pf — 2/3 and (b) free energy versus temperature for /i = 
0.84861. Note the large region of stability of the period-three 
phase in (a) and how the shape of the free-energy curves in 
(b) suggest that phase transitions to higher-period phases will 
be first order. 



as U increases further. The non-fermi-liquid character 
can also be seen in the self energy on the real axis. The 
imaginary part does have a quadratic behavior, but the 
curvature is the wrong sign and the intercept at u> = 
does not go to zero as T — > as it must in a fermi liquid. 
Similarly, the real part of the self energy is linear, but 
the slope has the wrong sign near to = 0. 

It is known in a fermi liquid with a local self en- 
ergy that the fermi surface and the DOS at the fermi 
level ( at T = 0) are unchanged as the correlations in- 
crease ilMuller-Hartmanrl Il989b|) . But because the FK 
mode l is not a fermi liquid there is no restriction on the 
DOS. Ivan Dongen and Leinund l)l997[) studied the Mott- 
like metal-insulator transition on the Bethe lattice in de- 
tail; in this study, one continues the homogeneous phase 
down to T — 0, ignoring any possible CDW phases. They 
found a simple cubic equation for the interacting DOS 
on the Bethe lattice which produces the results shown in 
Fig. One can see that for weak coupling, the DOS is 
essentially unchanged from the noninteracting case, but 
as U is increased further, the DOS first develops a pseu- 
dogap, where spectral weight is depleted near the chemi- 
cal potential, and then it gets fully reduced to zero. As U 
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when the correlation energy is large. We plot results for 
U = 2.5 where the localized electron DOS develops a 
pseudogap and then a gap region as T is lowered. Note 
how the gap develops as T is reduced below 0.4. 
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FIG. 12 Itinerant electron DOS for the spinless FK model on 
the Bethe lattice for different values of U ranging from 0.5 
to 3.0 in steps of 0.5 (p e = Pf — 0.5). The metal-insulator 
transition occurs at \J—2. It is preceded by a pseudogap phase 
for (1< U < 2). 



is increased further, the system splits into lower and up- 
per Hubbard bands, and the bandwidth of each subband 
decreases as U increases, while the separation increases 
between the two bands. This is precisely the kind of 
metal-insulator transition that Mott envisioned, where 
the DOS within a single band is suppressed to zero, but 
unlike Mott's prediction that the transition would gener- 
ically be a discontinuous first-order transition, here it is 
a continuous second-order transition. 

The metal-insulator transition occurs precisely at the 
point where the self energy develops a pole and diverges 
at uj = 0. This forces the DOS at the chemical potential 
to equal on any lattice. On lattices where the nonin- 
teracting DOS has band edges, the interacting DOS can 
remain zero within a correlation-induced gap. But on 
the hypercubic lattice, which has an infinite bandwidth, 
with an exponentially small DOS for large frequency, 
there is no precise gap, rather the DOS is exponentially 
small in the "gap region" and then becomes of order one 
within the Hubbard subbands. Indeed, on the hyper- 
cubic lattice, the DOS looks much like that in Fig. ^1 
but the metal- insulator transition occurs at U ~ 1.5 for 
Pe = Pf = 0.5. 

The DOS of the FK model does have a curious 
property — it is in dependent o f temp erature in the ho- 
mogeneous phase ijvan Dongenl ^992) (which was proved 
by mapping the FK model onto a coordination-three 
noninteracting Bethe lattice problem, which then has 
a temperature-independent DOS). The temperature- 
independence holds only in canonical formulations, where 
both p e and pf are separately fixed as functions of T [ig- 
noring the trivial shift of the DOS with p(T)]. In cases 
where the total electron concentration is fixed, but can 
vary between itinerant and localized electrons, the DOS 
will vary with T since pf generally varies with T. 

One can also e xamine the spectral function of the lo- 
calized electrons (jBrandt and Urbanekl Il992|h which is 
temperature dependent! The DOS shares similar be- 
havior with the conduction-electron DOS, having a gap 
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FIG. 13 Single particle DOS of the /-electrons for the spin- 
less FK model at half filling on a Bethe lattice with U = 2.5. 
Note how the gap opens up as T is reduced even though 
the calculations are above the CDW transition temperature 
(which occurs at T = 0.113746). The DOS is symmetric due 
to particle-hole symmetry. 



C. Falicov-Kimball-like Metal-Insulator Transitions 

The or iginal work of F al icov, Kimball, and 
Ram irez ijFalicov and Kimball Il969t iRamirez et all 
119701) studied a different type of metal-insulator 
transition — one where the character of the electronic 
states was unchanged, but their statistical occupancy 
fluctuated with temperature or pressure. As occupancy 
shifted from the itinerant to localized bands, the system 
underwent a charge-transfer metal-insulator transition. 
One of the main points of interest of this model was that 
the presence of a Coulomb interaction U between the 
two types of electrons could make the transition become 
discontinuous (first-order). The original approximations 
used a mean-field t heory which sh owed these first-order 
transitions. Later, iPlischkel l)l972j) proposed that when 
the coherent-potential approximation was applied to the 
FK model, the first-order transitions disappea r ed. T his 
was refuted bv lGoncalves da Silva and Falicovl l)l972j) . 

In the Falicov and Kimball approach to the metal- 
insulator transition, one works with spin-one-half elec- 
trons and fixes the total electron concentration (usually 
at ptotai — Pe + Pf — !)■ Alternatively, one can perform 
a partial particle-hole transformation on the /-electrons, 
and have a system where p e = pf and U is negative (at- 
tractive). There is a gap A to creating an electron-hole 
pair. We will concentrate on the electronic picture here 
(where one has Ef = U — A - 2t*) and the choice A = t* 
which most closely represents the original FK problem. 
We work on an infinite-coordination Bethe lattice. 

This model for metal- insulator transitions was so lved 
in infinite dimensions bv lChung and Freericksl l)l998|h It 
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FIG. 14 Conduction electron density as a lunction of 1/T for 
the spin-one-half FK model on a Bethe lattice with A = 1. 
The different curves correspond to different values of U. 



requires a fine tuning of the Coulomb interaction in order 
to have a first-order phase transition. This is illustrated 
in Fig. El where we plot the conduction electron density 
p e versus inverse temperature 1 /T for different values of 
U (recall Et also varies with U, since A is fixed at 1). 
The states that have p e — > as T — ► are insulators, 
while those with finite p e are metallic. We say a metal- 
insulator transition occurs when there is a discontinuous 
jump in the electron density, which goes from the metallic 
phase at high temperature to the insulating phase at low 
T. This occurs for 2.12 < U < 2.155 when A = 1. The 
shape of these curves is remarkably similar to those found 
in the mean-field-theory solution of lRamirez et all l)l97Cl() 
and verifies that first-order metal-insulator transition do 
occur within the CPA. 

The interacting electro nic DOS is plotted in F ig,IT5lfor 
four different values of U ijChung and Freerickslll99^ . In 
panel (a) we show the metallic case U — 2.16, where the 
DOS is large at the chemical potential (u> — 0) and does 
not depend strongly on T. As U is reduced to 2.15, where 
there is a discontinuous metal-insulator phase transition, 
we see a huge reconstruction of the DOS at the transi- 
tion, where the DOS at the chemical potential becomes 
small and decreases to zero as T — > 0. In panel (c), we 
see the evolution of the insulating phase for the strongly 
correlated insulator U = 2.12. Here the conduction elec- 
tron occupation does not evolve in a simple exponential 
fashion, but has a much sharper increase as T increases. 
The DOS has upper and lower Hubbard bands, with the 
lower band losing spectral weight as T — > 0. Finally, 
in panel (d), we show the DOS for a weakly correlated 
insulator. Here the conduction electron density evolves 
in an exponential fashion and there is no DOS at the 
chemical potential, since all of the electronic occupation 
is thermally activated. 



D. Intermediate-Valence 

The FK model does not have any hybridization be- 
tween the itinerant and localized electrons, so it can- 
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FIG. 15 Conduction-electron DOS for the spin-one-half FK 
model on a Bethe lattice with A = 1. The different panels 
are for different values of U, and the results are plotted for a 
number of different temperatures [the dashed lines are used 
in panel (b) for clarity] . The DOS has T dependence because 
pf varies with T. 



not have any quantum-mechanical mixed valence at finite 
temperature ijSubrahmanvam and Barmalll988^ which is 
quantified by having a nonzero average (c'f). Neverthe- 
less, one can have a classical intermediate-valence state, 
where the average /-electron occupancy lies between 
and 1 for T — > 0. In many such cases, one would expect 
CDW order (or phase separation) to take over in the 
ground state, but there are regions of parameter space 
where the system appears to remain in a homogeneous 
classical intermediate- valence state all the way to T = 0. 
We study such systems here. 

The intermediate-valence problem, and the possi- 
bility of the FK model having an instability to a 
spontaneo usly generated hybridizat ion, was first pro- 
posed by iPortengen et al\ (Il996albf) . Calculations in 
infinite-dimensions llCzvc holl. 1999) and one dimen- 
sion ijFarkasovskvl Il997l [1999) showed that a sponta- 
neously generated hybridization was unlikely over a wide 
range of parameter space, but they did not rule out the 
possibility everywhere. 

The first step to examining the intermediate-valence 
problem is to study the phase diagram of the FK model 
to see what different types of phases occur. We pick the 
spin-one-half model on the Bethe lattice again. We also 
pick p e + pf — 1 for the total electron concentration. 
Since the conduction electron fermi level is at when 
there are no /-electrons, if we pick —2 < Ef < 0, then 
the fermi level is pinned to Ef at U — 0, and one of two 
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FIG. 16 Intermediate- valence phase diagram for the spin-one- 
half FK model on a Bethe lattice with p e + pf = 1. Note 
how the classical intermediate valence state occupies only a 
small part of the phase diagram, being taken over by phase- 
separated, metallic, and insulating phases over most of the 
phase diagram. 



FIG. 17 Inverse of the hybridization susceptibility for the 
spinless FK model on a Bethe lattice with p e + pf = 0.5 
and Ef — —0.75. Note how the susceptibility monotonically 
increases with increasing U which suggests that it continues 
to diverge at T = 0. 



things can occur as U is increased: (i) either the classi- 
cal intermediate-valence state survives, or (ii) it phase 
separates into a mixture of a state with pf — and 
Pf = 1. The former occurs for Ef > —1 and the lat- 
ter for Et < — 1. As the interaction strength is increased 
further, then we find only two different stable solutions, 
a metallic phase, where pf — and an insulating phase 
where p* = 1. The zer o-temperature phase diagra m is 
summarized in Fig. HHI ifctiung and Freericksl l2000|) . In 
the phase-separated region, there are more possibilities: 
either the system remains phase separated in mixtures of 
two integer-valent states, or the phase separation has at 
least one intermediate-valence state in its mixture. The 
line of circles represents an approximate crossover line 
between these two possibilities (the intermediate-valence 
mixtures lie above the circles). Note how the metallic 
and insulator phases take over the phase diagram as U 
increases. In fact, the classical intermediate- valence state 
only occupies a small region of phase space, because 
it is unstable to phase separation over a wide region. 
This could explain why there are not too many classical 
intermediate- valence states seen in real materials. 

We conclude this section with a discussi on about the 
possi bility for spontaneous hybridization l)Zlatic et all 
2001). We confine our discussion to the spinless model 
with p e + pf = 0.5 for simplicity. If we repeat the above 
analysis, then the region of stability for intermediate va- 
lence at small U is the same as above —1 < Ef < 0. 
But as U is increased, it appears that the intermediate- 
valence state should be stable for all U when Ef — —0.75 
(but we have not performed a free-energy analysis to rule 
out the possibility of phase separation). We plot the in- 
verse of the hybridization susceptibility from Eq. I|135|) 
in Fig. El Note that it remains finite fo r all T. But 
we ca n get an analytic form when U — l)Zlatic et al\ . 
1200 1|) , which says that the susceptibility has a term pro- 
portional to |lnT|, which does diverge at T = 0. But 
this divergence does not guarantee that the ground state 



has spontaneous hybridization, in fact, it is known in this 
case that the ground state has (c^f) = ijFarkasovskvl 
2002), even though the susceptibility diverges. Here we 
find, down to the lowest temperature that we can calcu- 
late, that the hybridization susceptibility increases with 
U , which would imply that it continues to diverge at 
T = for all U. In this case, we do not know whether 
the ground state would possess a spontaneous h ybridiza- 
tion, or not, but the results of ISi et all l|l992f) suggest 
that spontaneous hybridization is indeed possible. In any 
case, since a real material will always have a nonzero hy- 
bridization (although it may be quite small), one could 
expect a quantum-mechanical intermediate valence tran- 
sition to occur at low temperature in regions of param- 
eter space where the hybridization susceptibility is large 
for the FK model. Alternatively, spontaneous hybridiza- 
tion may also occ ur if the /-electrons are allowed to 
hop l)Batistal. 12002]) . 

E. Transport Properties 

The optical conductivity, calculated from Eq. I|89|). is 
plotted in Fig. ^| for the spinless FK model on the hy- 
percubic lattice with w\ = 0.5. Panel (a) is the half-filled 
case (p e = 0.5) and panel (b) shows the behavior for dop- 
ing away from half filling. Both calculations are at low 
temperature T — 0.005 in the homogeneous phase (ig- 
noring any possible CDW phases). 

At half filling, the data behaves as expected. For small 
U, the system has a Drude like peak, whose width is de- 
termined by the scattering rate at low temperature (note 
the scattering rate does not vanish even at T = because 
the FK model is not a fermi liquid). As U increases, we 
see a charge-transfer peak develop, centered at v w U, 
and the low-energy spectral weight is suppressed because 
the system becomes a correlated (Mott-like) insulator 
(the metal- insulator transition occurs at U « 1.5 here). 
One can see that the transition is continuous, with Udc 
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FIG. 18 Optical conductivity of the spinless FK model with 
w\ = 0.5 on a hypercubic lattice at T = 0.005: panel (a) is 
the half filling case (p e = 0.5) with different values of U and 
panel (b) is the (7 = 2 case doped away from half filling. The 
numbers in panel (a) label the value of U, while the numbers 
in panel (b) show p e . 



smoothly approaching zero as the correlations increase. 

We concentrate on the U = 2 case as we dope away 
from half filling. One can see that as the system is doped, 
it becomes metallic, because the chemical potential now 
lies within the lower Hubbard band. As a result there 
is a transfer of spectral weight from the charge-transfer 
peak down to a Drude-like feature as the system is doped. 
What is remarkable, is that there is an isosbestic point 
present (which is defined to be a point where the optical 
conductivity is independent of doping and all the curves 
cross and typically occurs at v w U/2). The occurrence 
of such isosbestic points is present in a wide variety of 
models, but its origin is not well understood. 

We also plot the dc-conductivity, thermopower, and 
electronic contribution to the thermal conductivity, as 
derived in Eqs. (|93H98|I for a correlated system ((7 = 2, 
soli d line) with p K = 1 — w-\ an d five different w\ val- 
ues ijFreericks and Zlatio. l200 2b) . For these parameters, 
there is always a region of exponentially small DOS near 
the chemical potential at low temperature [but in this re- 
gion t{lo) decreases only as a power law] . As the localized 
electron concentration w\ moves away from 0.5, the high- 
temperature thermopower increases in magnitude due to 
the asymmetry in the DOS (it must vanish at 0.5 due to 
particle-hole symmetry) and the low-temperature ther- 
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FIG. 19 (a) DC conductivity, (b) thermopower, and (c) elec- 
tronic contribution to the thermal conductivity, for the spin- 
less FK model with p e = 1 — w\ and U = 2. Five fillings are 
shown: (i) W\ = 0.5 (solid line); (ii) wi = 0.4 (dashed line); 
(iii) wi = 0.3 (chain-dotted); (iv) u>i = 0.2 (dotted); and (v) 
w\ =0.1 (chain-triple-dotted). 



mopower shows a sharp peak for fillings close to half fill- 
ing (the sign is hole-like, because the DOS from the lower 
Hubbard band dominates the transport coefficients at 
low temperature); the dc conductivity and thermal con- 
ductivity both vanish at low T due to the "gap" as well. 
The thermoelectric figure-of-merit ZT — TajcS 2 / K e i is 
plotted in Fig. 120] — we find it is larger than one at high 
T for wi < 0.22, and for fillings close to half filling, there 
is a low-temperature spike in ZT that can become larger 
than one over a narrow temperature range. The spike at 
low T is due to the large peak in S and the small thermal 
conductivity; but the phonon contribution to the ther- 
mal conductivity can sharply reduce ZT if the phonon 
thermal conductivity is much larger than the electronic 
thermal conductivity (this all electronic calculation pro- 
vides only an upper bound to ZT). The Lorenz number 
is also plotted in Fig. [20| It gets huge at half filling, but 
becomes more metallic (« 7r 2 /3) as the filling moves fur- 
ther away from half filling. It is not a constant even at 
low temperature because the system is not a fermi liquid. 



F. Magnetic- Field Effects 

The magnetic field brings new features to the solu- 
tion, which we illustrate for the model with a fixed total 
number of particles and relatively large Falicov-Kimball 
interaction on a hypercubic lattice. We restrict the num- 
ber of /-particles per site to less than one, taking the 
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FIG. 20 (a) Lorenz number £k B /e 2 — K, e i/ad c T and (b) elec- 
tronic thermoelectric figure of merit ZT = Ta dc S 2 /K el for the 
spinless FK model with p e = 1 — wi and U = 2. Five fillings 
are shown: (i) Wi = 0.5 (solid line); (ii) w\ — 0.4 (dashed 
fine); (iii) wi = 0.3 (chain-dotted); (iv) wi = 0.2 (dotted); 
and (v) wi =0.1 (chain-triple-dotted). 



limit U** — > oo, and choose the position of the /- 
level such that there is a rather sharp crossover from 
the high-temperature state with a large concentration of 
/-electrons and a gap in the single particle DOS to the 
low-temperature state with a metallic conduction band 
(i.e., a fermi gas) and no /-electrons. That is, model 
parameters are such that the renormalized /-electron 
level is slightly above the chemical potential at T = 0. 
This is the regime that yields an an omalous magnetic re- 
sponse ijFreericks and Zlaticl Il998|) and is closely related 
to experimental materials like YbInCu4 which exhibit a 
valence-change transition. 

The average of the z-component of the /-electron mag- 
netization is 



Tr 



(137) 



where only the states with one /-electron in the presence 
of the magnetic field are considered. This lattice trace 
can be evaluated by using the cavity method again, where 
the lattice trace becomes equal to an impurity trace in 
the presence of an additional time-dependent dynamical 
mean field 



-'cavity 

z L 



Tr cf (e~ m ^ -» N -»f N ^S{\)m z f ) (138) 



where 5(A) is defined in Eq. (|104|l . The trace is per- 
formed by using the basis set which diagonalizes si- 
multaneously TL 



elements in the subspace without the /-particles and 
Zimp = Z L /Z cavity , we obtain 
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(139) 



where Zf = ^^^^[—^{Ef^H) — p,f}] is the parti- 
tion function of an isolated /-ion, Ef n (H) are the field- 
dependent eigenstates, and pj — p is the common chem- 
ical potential for the itinerant and localized electrons. 
The average magnetization is the product of a single-ion 
response and the average /-filling < pj < 1, i.e., the 
magnetic re s ponse of in dependent /-ions (Dzero, 2002; 
iDzero et all l2000t 12002(1 is modified by the interaction 
between the /-electrons and the conduction sea. This 
reduces the high-temperature Curie constant and intro- 
duces additional temperature and field dependences. 

We illustrate these features for the spin-1/2 model 
(s = 5 = 1/2) and plot in Fig.[2]thc temperature depen- 
dence of pf as a function of various parameters (includ- 
ing the magnetic field). The results can be explained 
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FIG. 21 The effect of the magnetic field on the /-electron 
concentration is shown for the spin 1/2 FK model with 
Ptotai = 1.5, U = 3, and Ef = —0.5, with gfiBH — 0.07 
(dotted line) and gp.sH = 0.17 (chain-dotted line). 

by noting that the renormalized /-level is just above 
the chemical potential and that a finite temperature in- 
duces an entropy driven "transition" (or crossover) into 
the magnetically degenerate state (the crossover temper- 
ature is denoted T v ). The field pushes the renormalized 
/-level closer to the chemical potential, which enhances 
the /-occupation, reduces the cross-over temperature and 
makes the transition sharper. For large enough fields, 
the concentration of /-electrons remains finite down to 
T=0, i.e., the system goes through a field- induced meta- 
magnetic transition. This is shown in Fig. 1221 where the 
magnetization obtained from Eq. i|139|) is plotted ver- 
sus magnetic field for various temperatures. Below the 
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FIG. 22 Localized electron magnetization for the spin 1/2 
FK model with ptotai = 1.5, U — 4 and Ef = —0.5 plotted as 
a function of magnetic field for various values of temperature, 
as indicated in the figure. 



crossover temperature, the low-field response is negligi- 
bly small (the Pauli susceptibility of the conduction elec- 
trons is neglected) but at large enough fields the mag- 
netization curves go through an inflection point, which 
indicates a crossover to a magnetic state (and is a meta- 
magnetic transition). Above the crossover, the curvature 
of the magnetization is positive and typical of a well- 
defined local moment. The effect of the magnetic field 
on transport properties is equally drastic. The magne- 
toresistance is plotted as a function of field in Fig. 1231 
and it shows that the low-temperature metallic state is 
destroyed above some critical field H c . Thus, the meta- 
magnetic transition in the /-subsystem is accompanied 
by metal-insulator transition in the conduction band. 
Taking the inflection point of the magnetization or the 
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FIG. 23 Magnetoresistance as a function of magnetic field 
plotted for several values of temperature. 

magnetoresistance, calculated for several values of U and 
Ef as an estimate of H C (T), we obtain the phase bound- 
ary which is shown in Fig. [21 together with the ana- 
lytic form H c (T)/H c (0) = y/l - (T/T*) 2 . The crossover 
temperature is renormalized by T* — T v /2. Note that 



the T* values in Fig. |^ differ by more than an order of 
magnitude, while the ratio kgT* / /isi7 c (0) is only weakly 
parameter dependent. 
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FIG. 24 Normalized critical field plotted as a function of re- 
duced temperature T/T* for several values of Ef and U. The 
solid line represents the analytic form \/l — (T/T*) 2 with 
T* = Tv/2. 



G. Static Holstein Model 



The stati c Holstein model l)Holsteinlll959tlMillis et all 
1995, 1996) can be viewed as a generalization of the FK 
model to the continuous-spin case. Like the FK model, it 
also displays metal-insulator transitions and CDW order 
phase transitions. We discuss both possibilities here. 

The phonon distribution function w(x) in Eq. I|39|l be- 
comes sharply peaked as T — > 0. If g ep is small enough, it 
approaches S(x), and the ground state is a noninteract- 
ing fermi gas (at half filling, the critical g ep is 0.8432 on 
the hypercubic lattice and 1.0854 on the Bethe lattice). 
Beyond this critical value of g ep , the phonon distribution 
function shows a double-peaked structure (becoming two 
delta functions at T = 0), and the "quasiparticles" scat- 
ter off of the local phonon even at T = 0. This creates 
a non-fermi-liquid state, and the DOS develops a pseu- 
dogap. As g ep is increased further, the pseudogap phase 
becomes fully gapped (when the self energy develops a 
pole), and a metal- insulator transition takes place. These 
results are illustrated for the spin-one-half model on the 
h ypercubic lat t ice in Fig. 1251 [similar to what was done 
in lMillis et al\ 111990) ] . 

The other area of interest is the CDW ordered phase 
at half filling. Since a static model does not supercon- 
duct, we need not worry about that order at all. One 
interesting puzzle in real materials is that the ratio of 
twice the CDW gap to the transition temperature is 
surprisingly large, usually much larger th an the BCS 
prediction of 3.5 ijBlawid and Millisl l200l() . The CDW 
transition temperature can be calculated by perform- 
ing ordered-phase calculations and extrapolating them 
to the point where the order d isappears to produce 
T c l|Ciuchi and de Pasaualel Il999j) . These calculations 
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FIG. 25 Plot of the electronic DOS (per spin) for the har- 
monic static Holstein model at half filling on the hypercubic 
lattice at T — 0. Different values of g ep have different thick- 
nesses of the lines. Note how the metal-insulator transition is 
continuous here too. 



show that even in the weak coupling limit, one finds a 
large ratio of 2A/fcsT c . Since CDW systems involve a 
distortion of the lattice to produce the ordered phase, 
one expects that anharmonic terms may play an impor- 
tant role in describing the physics behind them. So it is 
important to also examine what happens in the presence 
of anharmonic interactions. 

A surprising result was found when this system was 
analyzed on the hypercubic lattice. There was a uni- 
versal scaling law for the transition temperature, when 
T c was plotted against an extrapolated approximation to 
the wavefunction renormalization parameter (also called 
the quasiparticle Z factor, when one is in a fermi liquid), 
which is defined on the imaginary axis by 
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2 3^rT ' 



(140) 



The s caling law is plotted in Fig. 1261 (a) ((Freericks et all 
2000). The results for T c for a variety of coupling 
strengths and even for the case of quantum-mechanical 
phonons, all collapse onto the same curve. The results 
for an attractive Hubbard model (X symbols) do not, in- 
dicating that this scaling curve breaks down when the 
phonon frequency is made large enough. 

The other quantity of interest to examine is the ratio 
of the CDW order parameter at T — to the transition 
temperature. Since one might have expected the scaling 
theory for T c would have produced a universal plot there, 
that result might not have been so surprising. But, when 
anharmonicity is included, its effects should be stronger 
at low temperature, since the phonon distortion won't 
be allowed to be as large in the anharmonic case, and 
hence one would naively expect a smaller order parame- 
ter. But, as shown in panel (b), we see that the ratio of 
2A/ksT c still satisfies an a pproximate scaling law whe n 
anharmonicity is included ijFreericks and Zlaticl l2001a|) . 
Hence, one should expect, generically, that this ratio will 
be much larger than 3.5 except in extremely weak cou- 
pling cases (since it does approach 3.5 as g ep — > 0). 
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FIG. 26 Scaling plot of (a) the CDW transition temperature 
and (b) the CDW gap ratio versus the wavefunction renor- 
malization extrapolated at T c for the static Holstein model at 
half filling on the hypercubic lattice. In panel (a), we have 
also included some results for a quantum Holstein model, with 
a phonon frequency equal to 0.5, and for the attractive Hub- 
bard model (which is the infinite- phonon-frequency limit). 



IV. COMPARISON WITH EXPERIMENT 

A. Valence-Change Materials 

The FK model can be used to describe the anoma- 
lous features of rare earth intermetallic compounds that 
have an isostructural valence-change transition, as ob- 
served in the YbInCu4 and EuNi2(Sii_ a; Ge 2 ;)2 family 
of compounds. These intermetallics have been attract 



ing a lot of att e ntion recently llFelner and Novikl Il98f ; 
iFiguero^etaZjj 119981 ICarner et al\ |20 0(t iLevin et all . 



ll990HSarraoLll999HWada et aZ] . ll997tlzhang et oi.L l2002) 

and we describe the most typical features briefly. 

The temperature-dependent properties of YbInCu4, 
which have been studied most thoroughly and which we 
take as our example, are dominated at ambient pres- 
sure by a first-order valence-change transition at about 
40 K. The valence of Yb ions changes abruptly from 
Yb 3+ abo ve the transition temperature, T„, to Yb 2,85+ 
below T v l)Dallera et al 1 120021 iFelner and Novikl [l986) . 
The specific heat data shows at T v a first-order tran- 
sition with an entropy change of about AS ~ R In 8 
corresponding to a c omplete loss of magnetic degeneracy 
in the ground state llSarraol Il999). Neutron scattering 
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does no t provide any ev i dence for long-range order be- 
low T v ({Lawrence et all Il997|) . A t the transition, the 
/-occupation becomes non- integral l|Dallera et all l2002t) 
and the lattice expands by about 5%. The crystal struc- 
ture remains in the C15(b) class and the volume expan- 
sion estimated from the known atomic radii of Yb 3+ and 
Yb 2+ ions is compatible with t he valence change esti- 
mated from the Lrrr-e dge data ijCornelius et al\ . Il997t 
iFelner and Noviklll98^ . 

The low-temperature phase shows anomalies typical 
of a valence-fluctuating intermetallic compound. The 
electroni c specif i c heat and the susceptibility are en- 
hanced llSarraol 11999ft . the electrical resistance and 
the Hall constant are small and metallic, and the 
low-te mperature slope of the thermoelectric power is 
large llOcko and Sarraol l2002ft . The optical conduc- 
tivity is Drude like, with an additional structure in 
the mid-infrared r ange which appears qui te suddenly at 
T„ llGarner et all l2000|) . The ESR data ()Rettori et all 
1997) indicate a large density of states at the fcrmi level 
Ep . Neither the susceptibility, nor the resistivity, nor 
the Hall constant show any temperature dependence be- 
low T v , i.e., the system behaves as a fermi liquid with 
a characteristic energy scale Tfl 3> T v . The magnetic 
moment of the rare earth ions is quenched in the ground 
state by the f — d hybridization but the onset of the high- 
entropy phase cannot be explained by the usual Anderson 
model in which the low- and high-temperature scales are 
the same and the spin degeneracy is not expected to be 
recovered below Tfl- The valence-change systems, how- 
ever, recover the /-moment at T v which is a much lower 
temperature than Tfl- 

The high-T phase of YbInCu4 sets in at T v and is also 
anomalous. The Yb ions are in the stable 3+ configu- 
ration with one /-hole and with the magnetic moment 
close to the free ion value gp\/ J(J + 1)hb = 4.53/zb 
(c]l = 8/7 is the Lande factor and J = 7/2 is the total 
angular momentum of the 4/ 13 hole). The magnetic re- 
sponse is Curie- Weiss like with a small Curie- Weiss tem- 
perature 0, which does not seem to be connected with T„ , 
in any simple way l|Felner and NovikT , ll98(TtlSarracl , ll999fl 
(of course one does not expect a first-order transition 
temperature to be related to the Curie- Weiss tempera- 
ture, but this observation shows that the two phenomena 
are not governed by the same microscopic physics). The 
dynamica l susceptibility obtained from neu tron scatter- 
ing data l)Goremvchkin and Osbornl Il993|) is typical of 
isolated loc al moments, with w ell resolved crystal-field 
excitations ijMurani et all 12002') . However, neither the 
line shape of the dynamical response nor the tempera- 
ture dependence of the static susceptibility can be ex- 
plained by the Kondo model assuming Tk — T v . The 
Hall co nstant is large and n egative, typical of a semi- 
metal l)Figueroa et all Il998|) . the electrical resistance 
is also very larg e and not changed much by magnetic 
field up to 30 T l|lmmer et oi.L ll997'l. In typical Kondo 
systems, on the other hand, one expects a logarithmic 
behavior on the scale T/Tk and large negative mag- 



netore sistance. The high temperature optical conduc- 
tivity l)Oarner et all l2000l) shows a pronounced maxi- 
mum of the optical spectral weight at a charge-transfer 
peak near 1 eV and a strongly suppressed Drude peak. 
The high-temperature ESR data for Gd 3+ embedded in 
YbInCu4 resemble those fo und in integer-valence semi- 
metallic or insulating hosts ijAltshuler et oil Il995h . 

The anomalous properties of the high-temperature 
phase become most transparent if the valence-fluctuating 
phase is suppressed completely by pressure or dop- 
ing as in Ybi_ 2; Y 2 ;InCu4, where a substitution of 
15% Y ions s tabilizes the high-temperature phase 
down t o T=0 K llMitsuda et all l2002t lOcko and Sarraol 
l2002t IZhang et all l2002j) . The experimental re- 
sults for the resistivity, suscept ibility and the ther- 
mopower llOcko and Sarraol [2002) are shown in Fig. [Z7\ 
The susceptibility data show that the Curie- Weiss tem- 
perature for all the samples is about the same and much 
less than T v , i.e., the magnetic response of the high- 
temperature phase can be represented by a single uni- 
versal curve, provided one scales the data by an effective 
concentration of magnetic /-ions, p/, which is always 
smaller than the nominal concentration of /-ions. The 
functional form of the susceptibility above 10 K agrees 
well with the "single-ion" crystal field (CF) theory. At 
lower temperatures the susceptibility deviates apprecia- 
bly from the CF theory and shows a significant reduction 
of the Curie constant p/ which is not of the Kondo type. 

The resistivity in the high-temperature phase of 
Ybi_ x Y x InCu4 exhibits a weak maximum and the ther- 
mopower has a minimum above 100 K but neither quan- 
tity shows much structure at low temperatures, where the 
susceptibility drops below the single-ion CF values. That 
is, despite the presence of the well-defined local moments, 
there are no Kondo-like anomalies in transport or ther- 
modynamic properties of the high-temperature phase. 
On general grounds, one can argue that the disconti- 
nuity of the thermoelectric power at the valence tran- 
sition is a trivial consequence of the different thermo- 
electric properties of the two phases: the thermopower 
of the valence-fluctuating phase has an enhanced slope 
and grows rapidly up to T v , where it suddenly drops to 
the values characteristic of the high-temperature phase. 

The hydrostatic pressure, doping and the magnetic 
field also have, like the temperature, a strong effect 
on the properties of the valence-change materials. The 
critic al temperature of Y bInCu4 decreases with pres- 
sure ijlmmer et all Il997[) but the data are dif ficult to 
expla in with the Kondo volume collapse model l)Sarraol 
1999). Doping the Yb sites with Y 3+ or Lu 3+ ions 
reduces T v despite the fact that Y has a bigg er and 
Lu a smaller ionic radius l|Zhang et all l2002(k dop- 
ing the In sites by smaller Ag ions enhances T„ as 
in YbIni_. T Ag. 7 ,Cu4 fo r x < 0.3 l)Cornelius et all Il997t 
lLawrence et all |l999). Thus, doping cannot be ex- 
plained in terms of a chemical pressure. The low- 
temperature phase is easily destabilized by an external 
magnetic field: a critical field H C (T) induces a meta- 
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FIG. 27 Panel (a) shows the resistivity and the magnetic 
susceptibility of Ybi_ a; Y 2 :InCu4 as func tions of tempera- 
ture for various concentrations of Y ions (JOcko and SarraoL 
2002). Note, all the "high-temperature" data can be col- 
lapsed onto a single universal curve, by normalizing the sus- 
ceptibility with respect to an effective Yb-concentration (not 
shown). Panel (b) shows the thermopower of Ybi_ a; Y I InCu4 
as a function of temperature for various concentrations of Y 
ions (lOcko and SarraoL [2002ft . 



magnetic transition which is clearly seen in the magne- 
toresistance and the magnetization data. The experi- 
mental values of H C (T) define the H — T phase bound- 
ary and the analysis of many systems with different 
T° (zero-field transition temperature) and H c (0) (zero- 
temperature critical field) reveal the following universal 
behavior: H C (T) = H c (0)y/1 - (2T/T°) 2 . Despite the 
differences in T v and H c , the data give a constant ratio 
k B T°/fx B H c (0) = 1.8. 

In summary, the Yb systems switch at the valence- 
change transition from a low-entropy valence-fluctuating 
phase to a high-entropy magnetic semi-metallic phase. 
The transition is accompanied by the transfer of holes 
from the conduction band to the 4f shell and a metal- 
insulator transition (or crossover) . The high-temperature 
phase has degenerate local moments which interact with 
the conduction band (and, hence, get reduced) but the 
interaction is not of the usual Kondo type and no log- 
arithmic anomalies are seen. Such a behavior can be 
qualitatively understood by assuming that (i) the chemi- 
cal potential of the metallic phase is close to the location 
of the dip or gap in the DOS, which arises from many- 
body interactions, and that (ii) the energy of the localized 



magnetic configuration is above the non-magnetic one as 
T — > 0. The metallic phase is destabilized by magnetic 
fluctuations when the entropy gain due to additional (lo- 
calized) magnetic states overcomes the energy loss due 
to the change of the ionic configuration and the entropy 
loss due to the reduction of the number of holes in the 
conduction band. Once the chemical potential is brought 
close enough to the band edge, by doping, pressure, tem- 
perature or magnetic field, the entropy of the band states 
is suppressed and the magnetic entropy of the localized 
states becomes sufficient to compensate the energy loss 
and destabilize the metallic phase. This behavior is ob- 
served in many systems and it is not likely to be due to 
band-structure effects alone. Rather, the gap or pseudo- 
gap which is a necessary ingredient of the above scenario 
seems to be due to many-body interactions. 

The valence- change transition is also fo und in many 
Eu-b ased systems, and EuNi2(Sii_ x Ge x )2 ijWada et all 
Il997|) provides a typical and well-studied example. The 
thermodynamic anomalies are similar as in Yb systems 
but the transition temperature is higher and Eu ions un- 
dergo an almost complete valence change between the 
high-temperature, f 7 and the low-temperature f 6 con- 
figuration. That is, the transition is from a free-spin 
system at high temperatures to a simple metal at low 
temperatures. The electron transport is similar as in 
YbInCu4 but the intrinsic data are more difficult to 
measure because the volume change at the transition 
is large and samples typica lly crack when therm ally cy- 
cled. InEuNi 2 (Sii_ :c Ge a; )2 (|Sakurai et aZll200C|) the sign 
of the thermopower is reversed with respect to Yb sys- 
tems, which indicates an electron-like rather than hole- 
like transport. 

The high-temperature behavior of YbInCu4 and 
EuN^^ii-zGe^-like com pounds can be we l l de- 
scribed by the FK model (iFreericks and Zlaticl Il998t 
IZla,tic and FreeTk^. 1200114 12003 IZlatic et all 12001 jl. 
We describe the /-ions by two energetically close con- 
figurations which differ in their /-count by one, we take 
a common chemical potential for the conduction and lo- 
calized electrons (which is adjusted at each temperature 
to conserve the total number of electrons), we take the 
FK Coulomb repulsion between the conduction electrons 
(or holes) and the additional /-electron (or hole) large 
enough to open a gap in the conduction band. For realis- 
tic modeling we should also take into account the actual 
crystal-field structure including the splittings and the de- 
generacy of the ionic energy levels. For Yb ions, the state 
with no /-holes has unit degeneracy, while the single-hole 
case has a degeneracy of 8 (corresponding to J = 7/2). 
The two-hole case is forbidden due to the mutual repul- 
sion of two holes being too large. The 8-fold degeneracy 
of a single /-hole is further reduced by CF splittings and 
in a cubic environment we expect 4 doublets (unless there 
is some accidental degeneracy) . The external field further 
splits these CF states. For Eu ions, we take into consid- 
eration the 4/ 6 (3 + ) configuration with a non-magnetic 
ground state and two excited magnetic states, and the 
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magnetic 4/ 7 (2 + ) configurations. Since these states are 
pure spins states, there are Zeeman splittings but no CF 
splittings. All other states of Eu ions are higher in energy 
and neglected. 

To find the thermodynamic and transport properties 
we need the weights of the ionic 4/ configurations, i.e. we 
need the partition function of a more general FK model. 
Since the trace over / states still separates into a sum 
over the states which differ by one /-electron, the end 
result can be written as, 

Z imp = Z n o Y[ Z Qa (n) + Z n o ±1 TJ Z 0a (n - U) (141) 

cr a 

where ZQ a {ji) is defined by Eq. f I24JI and Z n o and Z^o^ 
are the partition functions of an isolated /-ion with 
and n°f ± 1 /-electrons, respectively. In the case of Yb 
ions, the — 14 configuration is non-magnetic, such 
that i? n o =1 4 = 1, and the rA — 1 configuration has one 

magnetic hole in the J = 7/2 spin-orbit state, such that 
the partition function becomes, 



77=1 



(142) 



The excitation energies Ef v (H) of a single magnetic hole 
are split into multiplets belonging to different irreducible 
representations of the crystal. 

In the case of Eu ions, the / 6 configuration has a non- 
magnetic ground state and two magnetic S — 1 and S = 
2 excitations. Assuming the usual Zeeman coupling with 
the magnetic field gives the ionic partition function, 
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FIG. 28 (a) Localized electron filling and (b) normalized 
spin susceptibility for eightfold degenerate (S = 7/2) Falicov- 
Kimball model with p t otai = 1-5, U — 2, and various Ef. 
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The / 7 configuration is the magnetic Hund's rule state 
(5 = 7/2) and its partition function is simply 

7/2 

Z„ 0=7 = e - pEs = 7 / 2 e-^ BHm '-^. (144) 

m a = -7/2 

The weights of the two ionic configurations are given by 
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which enter the Green's function formalism. 



(145) 



(146) 



In what follows, we discuss a case that is similar to 
that of Yb: the case of an S = 7/2 magnetic particle in 
the high-temperature phase and of a magnetically inert 
unoccupied state in the low-temperature phase; we ne- 
glect CF splitting but do consider Zeeman splitting. The 
generalization to Yb ions with the magnetic / 13 states 
split by the cubic CF, or to Eu ions with excited mag- 
netic states in the low-temperature / 6 configuration, is 
straightforward. We choose the Coulomb repulsion to be 
U = 2, which is large enough to open a small gap in the 
interaction DOS (see below) and we tune the bare /-level 
to bring the renormalized /-level slightly above the chem- 
ical potential. Thus, at T = 0, there are no /-electrons 
in the ground state but as T increases, the /-occupancy 
increases, producing a local-moment response. Since the 
filling of the / electrons is entropically driven, the filling 
increases rapidly with T as the degeneracy increases. In 
Fig. [2H1 we plot the average /-occupation and the spin 
susceptibility (normalized by the total m 2 z which equals 
42 for spin-seven-halves). Panel (a) of Fig. l2"51 shows that 
as Ef is reduced, the high temperature value of pf in- 
creases, the "transition temperature" decreases, and the 
transition becomes sharper. Panel (b) shows the spin sus- 
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FIG. 29 (a) DC resistivity and (b) thermopower for the (S = 
7/2) Falicov-Kimball model with ptotal = 1.5 and U = 2. The 
different curves correspond to different values of Ef. 
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FIG. 30 Interacting DOS for the (S = 7/2) Falicov-Kimball 
model with ptotal = 1-5, U = 2, £7/ = —0.6 and for various 
temperatures. 



occupation and the Falicov-Kimball interaction. 

The theoretical results shown in Figs. and [2D ob- 
tained for ptotal = 1-5, U = 2 and Ef close to —0.6, 
exhibit most of the qualitative features seen in the ex- 
perimental data. However, the model parameters and 
the absolute value of the temperature scale used for the 
static response functions cannot be determined unless 
one compares the dynamic properties of the model to 
experimental data. 



ceptibility, which is Curie-like but with a temperature- 
dependent concentration of local /-moments. Hence it 
has a peaked form, with a sharp reduction of the mag- 
netic response as T — > 0. Defining T v as the temperature 
at which the spin susceptibility drops to half of the max- 
imum value, we find T v = 0.05, 0.03, 0.007 and 0.001 for 
Ef = —0.4, —0.5, —0.6, and —0.65, respectively. In all 
these cases we are dealing with a crossover rather than 
a sharp phase transition. In the case Ef = —0.2 the 
peak is too broad for the transition to be defined, and 
for Ef = —0.7 the occupation of /-states remains finite 
down to lowest temperatures. A first-order phase transi- 
tion is possible for Ef between —0.65 and —0.7 but details 
of this transition cannot be obtained by the iterative nu- 
merical procedure, which becomes unstable. The effect 
of the /-level position on transport properties is shown in 
Fig. Panel (a) shows that the resistivity has a maxi- 
mum at an effective temperature T* which is much larger 
than T v and that the peak sharpens and moves to some- 
what lower temperatures as Ef is reduced. Effects at T v 
are not visible, except for Ef = —0.65, which lies close 
to the first-order transition. The thermopower results 
shown in panel (b) are somewhat similar, except that the 
transition at T v cannot be detected. In addition, there 
is a sign change in S that occurs at higher temperature. 
The transport anomalies are due to the development of 
a gap and the renormalization of the single-particle DOS 
at high temperature, driven by the increased /-electron 
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FIG. 31 Optical conductivity for the (S = 7/2) Falicov- 
Kimball model with ptotal = 1.5, U — 2, Ef — —0.6, and 
for various temperatures. 

The interacting conduction DOS is plotted in Fig.l30las 
a function of frequency, for the case ptotal = 1-5, U = 2, 
Ef = —0.6, and for various temperatures. Note, the 
temperature-induced shift in the chemical potential from 
the metallic region, with a high density of states, into the 
gap-region. Unfortunately, direct comparison with the 
photoemission spectra is not possible because the bulk 
and surface effects in valence-fluctuating systems are dif- 
ficult to separate. 
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Fig. shows the optical conductivity. We find a 
characteristic high-frequency hump at v = U, which 
is almost ^/-independent, and a temperature-induced 
transfer of spectral weight between the high- and low- 
frequency regions. A Drude peak grows as v — ► for 
T < T v . Thus, while the static transport is domi- 
nated by a pseudo-gap on the order of T*, the optical 
conductivity is sensitive to the valence-change transi- 
tion. The high-frequency peak in u(v) allows an esti- 
mate of the va lue of the FK Coulo mb interaction, yield- 
ing (7 « 1 eV (iGarner et all l2000j) for YbInCu 4 . Using 
T v = 42 K we find U/T v « 250 which can be used to 
narrow the choice of the model parameters. The experi- 
mental features can be reproduced by taking the /-level 
position at Ef = —0.6 but the procedure is not com- 
pletely unique, because U /T v is not independent of the 
total number of particles. As an additional constraint on 
the parameters one should demand that the theory repro- 
duces the peaks in the high-temperature transport prop- 
erties in the right temperature range. Clearly, a proper 
quantitative description should also incorporate the CF 
splittings of J=7/2 Yb states, which would provide the 
field-induced anisotropies of the magnetization. This has 
not yet been worked out in detail for either the Yb or Eu 
compounds. 



B. Electronic Raman Scattering 
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FIG. 32 Experimental Raman scattering re sults for (a) 
SmBe JNvhus et adll995blll997t) and (b) FeSi JNvhus et all 
1995a). The labels mark the different temperatures. 



Experimental Raman scattering data on a wide variety 
of different correlated insulators shows three distinctive 
features: (i) as the temperature increases, there is a sud- 
den transfer of spectral weight from a high-energy charge- 
transfer peak to lower energies; (ii) there is an isosbestic 
point, where the Raman response is independent of tem- 
perature at a special frequency, and all Raman scatter- 
ing curves cross; and (iii) if one takes the ratio of twice 



the spectral range where spectral weight is depleted at 
low temperature (representative of the insulating gap) to 
the temperature at which the low-energy spectral weight 
is restored (representative of T c ), then 2A/ksT c 3> 3.5. 
These features are shown in Fig.l32lw here we plot the Ra- 
man response for (a) SmBp iNvhus eiaH Il995bl fl997|) 
and (b) FeSi llNvhus et l!995alL 

Using the results of Tables fVTI and lyTTl at X = 1 allows 
us to calculate the Raman response in the A\ g and B± g 
channels. Since the Raman response ultimately is a com- 
plicated functional of the interacting density of states, 
and since all insulators have the same qualitative fea- 
ture of a gap near the chemical potential, we expect the 
Raman response to depend only weakly on the micro- 
scopic features of the insulating phase. Hence, we per- 
form our calculations for the simplest possible system 
that goes through a Mott-like metal-insulator transition: 
the spinless FK mod el at half filling (p e = w\ = 0.5) on 
a hypercubic lattice (IFreericks and Devereauxl l2001altJ : 
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FIG. 33 Theoretical calculation of the nonresonant Raman 
response in the (a) A\ & channel and the (b) Bi g channel at 
T = 0.5 and for various values of U . 

We begin with a plot of the Raman response at fixed 
temperature T — 0.5 and for a variety of values of U in 
the two symmetry channels fFig. l33f) . Note how for weak 
correlations we have a linear onset, and a higher-energy 
cutoff (from the band) as expected. As the correlations 
increase, we see the response separates into a charge- 
transfer peak centered at U and a low-energy peak (in 
the Bi g sector only). The vertex corrections in the Ai g 
sector suppress the low-energy features and we simply see 
the evolution of the charge-transfer peak, which is more 
asymmetric in shape. 

We choose to examine the inelastic light scattering for 
U = 2 in more detail, because it is a correlated insula- 
tor that lies close to the metal-insulator transition. We 
plot in Fig. the inelastic X-ray scattering response 
function along the zone diagona l for a variety of different 
temperatures ijDevereaux et qjl 12003^ . The response at 
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FIG. 34 Theoretical calculation of nonresonant inelastic X-ray scattering in the (a) A\ & channel and the (b) Bi g channel for 
U = 2 along the zone diagonal for a range of different temperatures. The thickness of the lines determines the temperature 
which ranges from 1.0 to 0.5 to 0.25 to 0.1. The individual curves for different X values have been shifted for clarity. The 
X = 1 curves correspond to conventional electronic Raman scattering. 



the zone center is the Raman response. We find that it 
displays all of the three features seen in experiment: (i) 
there is a low-temperature depletion of low-energy spec- 
tral weight, (ii) there is an isosbestic point, and (iii) the 
ratio of 2A/fcsT c (as extracted from the Raman response 
curves) is on the order of 10-20. These three features are 
generic to the plots at all X, except for X = 1 in the A lg 
sector, where the low-energy spectral weight disappears. 
This occurs because the finite-q vector mixes the differ- 
ent symmetries and the low-energy spectral weight is seen 
in the other symmetries. Note that the results are inde- 
pendent of symmetry channel for X — — 1, which makes 
for the possibility of an interesting experimental probe of 
nonlocal correlations. In general, we see a small amount 
of dispersion of the peaks and a generic broadening of 
the peaks as we move from zone center to zone bound- 
ary, but the isosbestic behavior remains for all q. So 
far resonant inelastic X-ray scattering experiments have 



been performed mainly at room temperature and with 
polarizers only on the incident beam of light, so direct 
comparison with theory is not yet feasible. 



C. Josephson Junctions 

Another area where the FK model has been applied 
to r eal materials is in the field of Josephson junctions 
[see iFreericks et al .] (l2002cft for a review]. A Josephson 
junction is a sandwich of two superconductors surround- 
ing a barrier material that can be a normal metal, a n 
insulator, or something in between ijJosephsonj. | l962). 
In ra pid single flux quantum (RSFQ) logic l|Likharevl 
2000), one tries to maximize the switching speed of the 
Josephson junction, while maintaining a nonhysteretic 
(single- valued) I —V characteristic. Since the integral of 
a voltage pulse over time is equal to a flux quantum for a 
Josephson junction, the height of the voltage pulse is in- 
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versely proportional to the width of the pulse; hence one 
wants to maximize the characteristic voltage to achieve 
the fastest switching speeds. The characteristic voltage 
is a product of the critical current at zero voltage I c with 
the normal-state resistance R n (slope of the I — V char- 
acteristic at high voltage). 

In conventional tunnel junctions, the barrier material 
is an insulator, so I c is low and R n is high; in proximity- 
effect junctions, the barrier material is a normal metal, 
so I c is high and R n is low. Is it possible that one can 
maximize the product of I c R n by choosing a material 
to lie close to the metal-insulator transition, where both 
I c and R n can be large? We can examine this question 
by describing the barrier material with the spin-one-half 
FK model on a cubic lattice in the local approximation 
with wi = 0.5 and p e = 1. This system has a metal- 
insulator transition in the bulk at U » 4.91 To sim- 
ulate the properties of a Josephson junction, we must 
solve an inhomogeneous DM FT problem as first done by 
iPotthoff and Noltind (l999) and generalized to th e su- 
perconducting state bv iMiller and Freericksl l|200llh Us- 
ing the FK model to describe the barrier material was 
performed in lFreericks et al\ ll200l"bll2003tl and has been 
reviewed elsewhere ()Freerickse7T^.U2002c|) . 




2 4 6 

Interaction Strength U [t*] 

FIG. 35 Figure of merit as a function of the FK correlation 
strength in a Josephson junction with a barrier of (a) one and 
(b) five planes. The dotted line is the Ambegaokar-Baratoff 
prediction for an ideal tunnel junction. The plots are at two 
temperatures T = T c /ll (open symbols) and T = T c /2 (filled 
symbols). A(0) is the superconducting gap at T = 0. 

In Fig. 1351 we plot the figure of merit I c R n versus U 
for two cases: (a) a single-plane barrier and (b) a five- 
plane barrier. Notice how the I c R n is maximized in a 
ballistic metal for iV = 1 at low T, but then is maxi- 
mized in the tunnel junctions as T increases. The flat- 
ness o f the curves in panel (a) for U > 5 is a verification 
of the lAmbegaokar and Baratofll l)l963f) analysis which 



says I c R n is independent of the properties of the barrier 
for thin tunnel junctions. Our result is somewhat lower 
than the Ambegaokar-Baratoff prediction due to fermi- 
surface effects, proximity and inverse proximity effects, 
etc. What is interesting is that in panel (b) there is a 
marked increase in the figure of merit at low temperature 
near the metal-insulator transition. Indeed, we find I c R n 
is maximized on the insulating side of the metal-insulator 
transition, and the optimization remains for a wide range 
of temperature (the reduction in I c R n arises mainly from 
the temperature dependence of R n for a correlated in- 
sulator). These results are consistent with experiments 
performed on junctions made out of N bTiN for the super - 
conductor and Ta^N for the barrier l|Kaul et all 12001). 
As tantalum is removed from TaN, it creates tantalum 
vacancies, which are strongly interacting with the con- 
duction electrons and can trap them at the vacancy site. 
Such physics is described well by the FK model, with the 
tantalum vacancy sites serving as the "/-electrons" that 
scatter the conduction electrons. 
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FIG. 36 Critical current as a function of temperature for a 
variety of different Josephson junctions. Note how the corre- 
lated insulator (solid triangles) is parallel to the thin insulator 
(solid line) for a wide range of T. 



In Fig. 1361 we plot the critical current as a function 
of temperature. Once again, our c alculations produce 
the lAmbegaokar and Baratofj l)l963|) results for thin tun- 
nel junctions. As junctions are made more metallic, or 
thicker, the critical current drops more rapidly as a func- 
tion of temperature. But the correlated barrier, just on 
the insulating side of the transition, has a thermal slope 
dI c (T)/dT that is essentially equal to that of the thin 
tunnel junction in the range 0.4T C < T < 0.7T C which is 
a typical operating range for a junction. Hence these cor- 
related metal barriers (or SCmS junctions) may be the 
optimal choice for the barrier in a Josephson junction, 
enabling ultrafast superconductor-based digital electron- 
ics. 
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D. Resistivity Saturation 

The semiclassical Boltzmann equation approach to 
transport indicates that the resistivity of a material 
should continue to rise as the temperature rises. How- 
ever, the theory breaks down once the mean free path of 
the electrons becomes shorter than the interatomic spac- 
ing. Hence, there is an expectation for the resistivity to 
slow its increase once the mean free path becomes too 
small. An interesting set of materials that have been in- 
vestigated a quarter century ago are t he so-called A15 
compounds such as VaSi and NbjGe ijFisk and Webbt 
1976). These materials sparked much interest as being 
the highest temperature superconductors of their day, 
and they were widely studied, but a number of features 
about these materials remain unsolved. 

In conventional electron-phonon scattering metals, the 
resistivity behaves like 

p(T)=A p T + B p (147) 

at high temperature. Here A p is proportional to the 
electron-phonon coupling strength and B p is proportional 
to the impurity concentration. What was found in the 
A15 materials was that A p was much smaller than ex- 
pected, given the known strength of the electron-phonon 
coupling and B p was sizable, even in very pure samples. 
This behavior was called resistivity saturation, since it 
implies a slow turnover of the high-temperature resistiv- 
ity 

I 1999. iMillis et all ||1999|) solved an important puz- 
zle in the resistivity saturation problem. They examined 
the static harmonic Holstein model at high temperature 
and saw the characteristic shape seen in resistivity sat- 
uration for strongly coupled systems. This analysis ne- 
glects the quantum- mechanical nature of the phonons, 
since their kinetic energy is neglected, but this should be 
a good approximation if the phonons are at a tempera- 
ture much higher than the Einstein frequency of the Hol- 
stein model (estimated to be about the Debye frequency 
of a real mater i al) . N ote that the resistivity formula used 
m Millis et al 1 (I1999D does not properly have the velocity 
factors for the Bethe lattice. This does not change the 
results much near half filling, but would have a larger 
effect for fillings closer to the band edge. 

We illustrate this phenomenon in Fig. [37| where we 
plot the resistivity at half filling p e = 1 for the spin-one- 
half static harmonic Holstein model on the hypercubic 
lattice. We examine coupling strengths ranging from a 
weakly coupled system (that is a fermi gas at T = 0) to 
a strongly coupled insulator, that develops a gap in the 
conduction electron DOS (see Fig. l2l>)l . 

The theory behind the resistivity form in Eq. I|147ll de- 
pends on three assumptions: (i) the temperature is much 
higher than the phonon energy scale, so the ions can be 
approximated by classical oscillators and the equiparti- 
tion theorem says (xf) ~ T/k; (ii) the electron-phonon 
interaction can be treated in second-order perturbation 
theory so the scattering rate is g1 p {xj); and (iii) Boltz- 




FIG. 37 (a) Resistivity as a function of temperature for the 
half filled static harmonic Holstein model on the hypercubic 
lattice and (b) parameters A p (solid line) and B p (dashed line) 
fit for the range of temperature 0.2 < T < 1. Note how there 
is a wide range of temperature where the resistivity appears 
linear, but with a large intercept. Note further how A p grows 
slower than gl p oc A. 



mann transport theory can relate the scattering rate 
to the resistivity. Assumptions (i) and (ii) imply that 
the sc attering rate is linear in temperature. IMillis et all 
( 1999) find that the scattering rate actually increases like 
VT due to the breakdown of the second-order perturba- 
tion theory for strong coupling. This is the key to re- 
sistivity saturation, as can be seen in Fig. 1371 Indeed, 
forcing the square root behavior to fit to a linear form 
necessarily produces all of the observed resistivity satura- 
tion phenomena. Note further that even for an insulating 
phase, we see the same kind of "resistivity saturation" at 
high enough temperatures. 

It is interesting as well that there are a variety of mate- 
rials (mainly the doped fullerenes and the high-T c super- 
conductors) where the resistivity obeys the linear form 
of Eq. (|147fl with small or vanishing B p up to very high 
temperatures, where a naive estimate of the mean free 
path is much less than a lattice spacing. The DMFT 
analysis given here does not shed light onto that puzzle, 
but rece nt work using other techniques ha s made much 
progress l|Calandra and Gunnarssonll200l[) . 
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E. Pressure-Induced Metal-Insulator Transitions 

Nickel iodide (N1I2) is a transition-metal halide 
that undergoes an isostructural metal-insulator tran- 
sition at room temperatu re as a function of pres- 
sure ({Pasternak et aZJ .fl990"l . Nickel iodide crystallizes in 
the CdCl2 structure, which consists of alternating hexag- 
onal planes of nickel and iodine. The nickel ions have 
a 2+ valence and the iodine ions have a 1— valence; 
the neutral sandwich I~-Ni ++ -I~ is stacked vertically to 
form the crystal. In the insulating state, the 5 = 1 
nickel spins order in a helical spin density wave, that 
is closely approximated by ferromagnetic nickel planes, 
stacked in an antiferromagnetic fashion. At ambient pres- 
sure Tjv = 75 K, and Tjy increases by a factor of four to 
310 K at 19 GPa, where the system undergoes a metal- 
insulator transition and becomes metallic. The physical 
picture is as follows: the d-electrons in the nickel band are 
strongly correlated and form a Mott-like insulating state 
(which can be approximated by dispersionless localized 
electrons), and the iodine p-bands are completely filled 
by the transfer of an electron to each iodine. As pres- 
sure increases, the relative position of the iodine p-bands 
and the nickel d-bands changes, with the iodine p-bands 
moving closer to the fermi level. Once the p-bands reach 
the fermi level, electrons spill from the p-bands to the 
nickel d-bands, quenching the magnetic moment (chang- 
ing Ni ++ to Ni + ), stopping the antiferromagnetic tran- 
sition, and allowing hole conduction within the p-bands. 
Hence there is a transfer of charge from the iodine to the 
nickel as the pressure increases causing a metal-insulator 
transition. This is precisely the kind of transition envi- 
sioned in the original Falicov-Ki mball model and o ptical 
experiments confirm this picture ijChen et alJ . ll993|h 

Indeed, one can add spin interaction terms to the FK 
model, to more accurately approximate this system, and 
perform an analysis of the thermodynamics and map 
out a phase diagram. This has a lready been done with 
a me an-field-theory calculation ijFreericks and Falicovl 
1992), which shows antiferromagnetic insulating phases, 
paramagnetic insulating phases, and metallic phases and 
is depicted in Fig. |2H1 The phase diagram also shows 
a classical critical point, where the first-order metal- 
insulator transition disappears. Unfortunately, the criti- 
cal point is estimated to lie around 1400 K, which is be- 
yond the disintegration temperature for N1I2. It would 
be interesting to repeat the analysis of MI2 with DMFT, 
to see if the estimate of the critical point is reduced in 
temperature and to examine properties like optical con- 
ductivity which have been measured. 



V. NEW DIRECTIONS 

A. 1/d Corrections 

The dynamical mean field theory is exact in the 
infinite-dimensional limit. But the real world is finite 
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FIG. 38 Mean-field theory phase diagram for Nilj 
calcul ated within the Falicov -Kimball model pic- 
ture flFreericks and Falicovl |l992). The shaded region 
is the antiferromagnetic insulator, the region above the 
chain-dotted line is the paramagnetic insulator. The param- 
agnetic metallic phase lies to the right of the nearly vertical 
solid line. The first order transition line from an insulator 
to a metal ends in a classical critical point (above that 
temperature, one can continuously crossover from the metal 
to the insulator). The symbol A denotes the difference 
in energy from the localized Ni electrons to the bottom 
of the I bands. It varies with an external parameter like 
pressure. The dotted lines denote the regions of phase space 
where three solutions of the mean-field theory equations are 
stabilized. 



dimensional, and so it is important to understand how 
the infinite-dimensional limit relates to finite dimensions. 
The simplest approximation for finite dimensions that 
one can take is the so-called local approximation, where 
one performs the DMFT, but uses the correct finite- 
dimensional noninteracting DOS in the relevant Hilbcrt 
transform . Such an approa ch was adopted for the FK 
model bv iFreericksl (|l993b). Since we expect Hartree- 
Fock theory to be asymptotically exact in the limit as 
U — > (for d > 1), and since the DMFT uses the correct 
noninteracting susceptibility in determining CDW order, 
the local approximation is asymptotically exact at weak 
coupling. Deviations are expected to be much stronger in 
the large coupling limit (where the model can be mapped 
onto an effective Ising model at half filling), due to the 
renormalization of T c from spatial spin- wave fluctuations 
that arc neglected in the local approximation. In fact, 
the situation at weak coupling is somewhat more com- 
plicated, because it is well known that the Hartree-Fock 
approximation to the Hubbard model is renormalized due 
to local qua ntum fluctuations by facto rs on the order of 
three to five (|van Dongenlll991blll994T) . The 1/d correc- 
tions to this renormalization are typically small for the 
Hubbard model (on the order of a few percent in two 
and three dimensions). This weak-coupling perturbative 
analysis has not been carried out for the FK model, and 
there may be surprises there, because the DMFT T c at 
weak coupling does not take the conventional exponen- 
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tial form. But the general thrust of this analysis is that 
one expects DMFT and the local approximation to be 
reasonably good at weak coupling, with larger deviations 
as the coupling is made stronger. 

The first attempt at a systematic expansion in X/d 
th at can include nonpe r turba tive effects was made 
by ISchiller and Ingersentl (0.995). They constructed a 
generic self-consistent two-site impurity problem, which 
had feedback into an auxiliary single-site impurity prob- 
lem to calculate both the local and nearest-neighbor con- 
tributions to the self energy (i.e. the self energy is al- 
lowed to have an additional momentum dependence pro- 
portional to ek). This approach suffered from two diffi- 
culties. First, their self-consistent equations encountered 
convergence problems when U was made large, and they 
were only able to achieve converged results for small U 
(where one does not expect there to be large corrections) . 
Second, one can immediately see that in regions where 
the local contribution to the self energy is small, one 
might lose causality (the imaginary part of the self en- 
ergy can have the wrong sign) if the coefficient of the 
ek term is too large. Indeed, the numerical calculations 
suffered from negative DOS at the band edges as well. 

A number of different approximate methods were also 
attempted to improve upon the situation. The first 
method is called the dynamical cluster approximation 
(DCA) pettier et all l2000l Il998j) and we will describe 
it in detail below. The second one is based on a trunca- 
tion of a memory function expansion (which uses the Li- 
ouville operator space and conti nued-fraction techniques 
to de termine Green's functions) ijMinh-Tienl Il998l iTranl 
1999). This technique is approximate and it includes 
static nonlocal correlations and dynamic local correla- 
tions, which results in some pathologic behavior, such as 
the interacting DOS remains temperature independent 
in the canonical ensemble, just like DMFT. The third 
method is based on a moment analysis of the Green's 
function and uses a self-consistent feedback of the /- 
electron charge susce ptibility onto the conduction elec - 
tron Green's function ijLaad and van den Bosschell2000|) . 
This technique also requires a number of uncontrolled ap- 
proximations in order to have numerical tractability. One 
of the controversial predictions of the moment approach 
is that the FK model in two dimensions has the Mott-like 
metal-insulator transition occur at U = 0. The first two 
techniques predict that it occurs at a finite value of U on 
the order of half of the bandwidth. 

We will concentrate our discussion here on the dynam- 
ical cluster approximation, which is a systematic tech- 
nique in the thermodynamic limit for incorporating non- 
local correlations into the many-body problem. It is not 
a formal expansion in 1/d, but rather is an expansion in 
1/N C , where N c is the size of the self-consistent cluster 
employed in the computational algorithm. From a phys- 
ical point of view, one should view the DCA as an ex- 
pansion in the spatial size over which spatial fluctuations 
are included, so it allows large momentum (short-range) 
spatial fluctuations, but does not properly describe long 



wavelength fluctuations. The basic idea is to allow the 
self energy to have momentum dependence in a coarse- 
grained fashion. The Brillouin zone is divided into N c 
equally sized regions, and the self energy assumes con- 
stant values within each of these regions (but can vary 
from region to region). It turns out that one can guar- 
antee causality is preserved with this technique. 

We now describe the DCA algorithm in detail. We be- 
gin by dividing the Brillouin zone into N c coarse-grained 
cells, labeled by the central point of the cell K. In this 
fashion, every wavevector in the Brillouin zone can be 
written as k = K + k with k ranging over the coarse- 
grained cell. The K points are chosen to correspond 
to the wavevectors of the A c -site cluster with periodic 
boundary conditions. The algorithm proceeds then as 
follows: (i) we choose our initial guess for the coarse- 
grained self energy S(K,o;); (ii) we construct the coarse- 
grained Green's function via 

G(K, W ) = ^]T— — ^— ; (148) 

N ^ uj + \i - S(K, uj) - e K+ic 

(iii) we extract the effective medium for the cluster 
Go^K,^) = G-^Kjw) + E(K,w); (iv) we solve the 
cluster problem for the cluster Green's functions G c given 
the effective medium Go and determine the cluster self 
energy via S c (K,cj) = Gq~ (K, w) — G f T 1 (K,cj); (v) we 
equate the cluster self energy with the coarse-grained self 
energy S(K,w) = S C (K, u>) and substitute into step (ii) 
to repeat the process. The algorithm is iterated until it 
converges. Once converged, one can then evaluate the 
irreducible vertex functions on the cluster using analo- 
gous Dyson equations; equating the cluster irreducible 
vertex functions to the coarse-grained irreducible vertex 
functions on the lattice then allows one to compute sus- 
ceptibilities. 

Under the assumption that the technique used to solve 
the cluster Green's functions given the effective medium 
Go is a causal procedure, then the cluster Green's func- 
tions and self energies are manifestly causal. Forming 
the coarse-grained Green's function on the lattice from 
the causal cluster self energy maintains the causality of 
the Green's function. The only place where noncausality 
could enter the algorithm is in the step where we ex- 
tract the effective medium, if the imaginary part of G _1 
is smaller in magnitude than the imaginary part of the 
coarse-grained self energy. But one can prove that this 
never occurs, hence the algorithm is manifestly causal. 

The DCA was applied to the FK model using two 
different techniques. When the cluster size was small 
enough (N c < 16), one could use exact enumeration to 
determine the weights for every possible configuration of 
localized electrons on the lattice, and since the action is 
quadratic in the fermions, one can determine the parti- 
tion function exactly. For larger clusters, one needs to 
determine the weights via a statistical sampling proce- 
dure, including the weights that are largest in the parti- 
tion sums and neglecting those that are too small. We 
will describe the exact enumeration method only here. 
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We let / denote a configuration n{ , n f N } of 
the localized electrons (we consider only the spinless FK 
model for simplicity). Then, the partition function for 
the cluster becomes 

{/} 

x -Q Dct[G ' 1 (z, j,iu n ) -Un{6jj] 



noj, 



where the symbol N[ denotes the total number of lo- 
calized electrons in the configuration / and i and j are 
spatial indices on the cluster. The determinant is over 
the spatial indices i and j. The weight for the configura- 
tion / is the corresponding term in Eq. i|149l) divided by 
Z c ; we denote that weight by w/. The cluster Green's 
function on the imaginary or real axis is then simply 



G c (i,j,z) = ^2w f [G 1 {i,j,z) - Un{Sij]' 
{/} 



(150) 



The weights are determined solely from the Green's func- 
tions evaluated on the imaginary axis, while the Green's 
functions on the real axis are trivial to determine once 
the weights are known. The convergence of the iterative 
algorithm is much slower on the real axis though. 
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FIG. 39 Many-body density of states for the half-filled FK 
model on a s quare lattice using t he DCA on a 4 x 4 cluster 
with U = 4 ijHettler et all 1200(f) . Note how the DOS now 
has temperature dependence, and how a pseudogap develops 
and deepens as T is lowered. Once a CDW opens, the system 
will open a true gap. 

One of the interesting features of the DCA is that 
it restores temperature dependence to the conduction- 
electron DOS. This is shown in Fig. 1391 where we plot the 
DOS at five different temperatures (all above the CDW 
T c ) for an intermediate value of U = At on a 4 x 4 square- 
lattice cluster at half filling ijHettler et all 120 00). Note 
how the pseudogap deepens as the temperature is low- 
ered. Of course a true gap will develop in the CDW 
phase at the lowest temperatures. There is no obvious 
way of "turning off" the CDW phase below T c anymore, 
because the / configurations with the highest weight will 
be those that correspond to the CDW ordered state be- 
low T r . 
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FIG. 40 Phase diagram for the half-filled FK m odel on a 
square lattice with the DCA jHettler et aZ.ll200CD . The dif- 
ferent lines correspond to different size clusters. The 2x2 
cluster determines T c from the ordered phase, while the other 
clusters employ a susceptibility analysis. Also included in the 
figure is the Ising model T c for the effective spin model with 
J = t 2 /2U. 



We also plot the phas e diagram in Fig. UJJat half fill- 
ing ijHettler et all 12000^ . Included in the plot are the 
local approximation, cluster calculations from suscepti- 
bilities on clusters up to N c = 36, and an Ising model T c 
for J = 2t 2 /U which the FK model maps onto at strong 
coupling; note that the 2x2 lattice result comes from an 
ordered-phase calculation with exact enumeration, rather 
than a susceptibility calculation. There are a number of 
interesting comments to make about this phase diagram. 
First, as expected we find the T c is hardly changed on 
the weak-coupling side of the diagram, indicating that 
the local approximation determines the majority of the 
many-body correlations in this regime. In the large-t/ 
regime, one can see that the DCA does converge towards 
the correct Ising model result, albeit slowly. The 2x2 
data is surprising, as it is much lower than the results 
for other clusters, but one can perhaps understand why 
this would be by realizing that a four-site square lat- 
tice cluster with periodic boundary conditions is identi- 
cal to a four-site one-dimensional chain cluster, except 
for a change in the overall normalization of the hopping. 
Since the T c for a one-dimensional chain is equal to zero, 
the suppression of T c for the four-site cluster could be 
arisi ng from the two pos sible interpretations of the clus- 
ter. iHettler et all l)2000|) argue in a similar vein, where 
they observe that the fermi-surface points on a four-site 
cluster occur at the van-Hove singularities, which have 
uncommonly large scattering rates and hence are unfa- 
vorable for CDW order. 

It is possible to explore dynamical susceptibilities 
within the DCA as well. One simply needs to extract the 
frequency and coarse-grained momentum irreducible sus- 
ceptibilities and use the relevant Bethe-Salpeter equation 
to calculate the full susceptibility. Unlike in the DMFT, 
where the irreducible vertex only renormalizes symmetry 
channels that have the full symmetry of the lattice, the 
vertex now renormalizes all channels, because it contains 
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momentum dependence. This remains one of the open 
problems of interest for the DCA. 



B. Hybridization and /-Electron Hopping 

The FK model in infinite dimensions admits a sim- 
ple solution, because the Hilbert space of the impurity 
problem generated by the DMFT mapping is a direct 
product of invariant subspaces with a fixed /-occupation. 
In each subspace we can easily find the itinerant elec- 
tron Green's function and the generalized partition func- 
tion [the generating functional of Kadanoff and Baym, 
Z = Dct(G " 1 )]. In the absence of quantum fluctuations 
in the /-particle number each /-state has an infinite life- 
time but the /-spectral function acquires a finite width 
due to the / — c coupling which acts as an effective time- 
dependent potential for the /-electrons and the broad- 
ening is due to the statistical averaging over all possible 
states. In this respect the FK problem is a lattice gen- 
eralization of the X-ray edge problem. Actually, as long 
as the coherent scattering of conduction electrons on /- 
electrons is neglected, we can drop the self-consistency 
condition Eq.(31) and use the X-ray e dge model solution 
for the lattice problem l)Si et a/lll992(l . At lower temper- 
atures, the system becomes coherent and the single-site 
X-ray solution has to be replaced by the DMFT self- 
consistency (effective medium) , so as to keep track of all 
the other /-sites on the lattice. Despite its simplicity, 
the DMFT solution of the FK model has some interest- 
ing features and can be used to describe physical systems 
in which quantum fluctuations can be neglected, like the 
high-temperature phase of Yb- and Eu-compounds with 
a valence- change transition or the charge-transfer metal- 
insulator transition in N1I2. 

However, if the /-particle dynamics are important, 
which is often the case at low enough temperature, one 
should consider a more general model. For example, in 
the Yb-systems discussed in Section IIV.AI the valence 
of Yb ions below T v is much larger than predicted by 
the FK model and the reduction of the local moment 
does not follow from Pf(T) but is due to quantum fluc- 
tuations and lifetime effects. If the ground state is a 
mixture of 2+ and 3+ states one should not neglect the 
f—c hybridization and a better description would be pro- 
vided by a periodic Anderson model with an additional 
FK term. The actual situation pertaining to Yb ions is 
quite complicated, as one should consider an extremely 
asymmetric limit of the Anderson model, in which there 
is no Kondo resonance (the ground state is not Kondo- 
like) and there is no single u niversal energy scale w hich is 
relevant at all temperatures l)Krishnamurtv et al\ .\l980). 
We believe such a generalized model would behave as 
the FK model at high temperatures and as a periodic 
Anderson model at low temperatures. Indeed, incoher- 
ent scattering of conduction electrons on the FK ions 
favors the gap opening even in the presence of hybridiza- 
tion; and if the width of the /-level is large due to the 



FK interaction, the additional effects due to quantum 
mixing should be i rrelevant. This is supporte d by the 
scaling arguments l|Withoff and Fradkinl 11990^ and the 
numerical reno rmalization group anal ysis of various im- 
purity models l|Voita and B ulla. 2002j) which show that, 
as long as the chemical potential is close to the gap and 
the hybridization is below some critical value, the /-level 
decouples from the conduction sea. On the other hand, 
at low temperature, the chemical potential shifts away 
from the gap in the DOS. For a flat conduction band 
and the /-level close to the chemical potential, the hy- 
bridization effects are important and they can drive the 
system towards a valence- fluctuating fixed point. The 
most likely effect of the FK correlation at low tempera- 
tures is to renormalize the parameters of the Anderson 
model. 

A generalized FK model might also be needed 
to describe the electron ic ferroelectrics proposed by 
iPortengen et all lll996albT ). The Bose- Einstein conden- 
sation (BEC) of / — c excitons and a spontaneous po- 
larization can be obtained by a mean-field treatment of 
the hybridized model. This picture is further supported 
by the examination of the spontaneous hybridization for 
the spinless FK model on the Bethe lattice, where it is 
found that as the temperature is lowered the system ap- 
pears to have a logarithmic divergence in the spontaneous 
hybridization susceptibility at T = 0. Normally we can- 
not reach such a state because the system will have a 
phase transition to either a phase separated state or a 
charge-density wave, but we can tune the system so that 
it remains in a classical intermediate- valence state down 
to T = 0. When this occurs, effects of even a small hy- 
bridization will take the system away from the FK fixed 
point at low enough temperature. 

Another generalization is obtained by allowing direct 
/ — / hopping, t{p which leads, for t{ 4 = t ? -, to t he rich 
physics of the Hubbard model l)Hubbardl.tl963j) . The 
case tL <C ty is also interesting, because even a small 
f — f hopping induces an / — c coherence and gives 
a ground state (for large U) which is either a BEC of 
electron-hol e pairs (/ — c excitons) or an orbitally or- 
dered state llBatistall2002l) . Unfortunately, as soon as the 
/-electron dynamics is restored, the evolution operator 
connects all sectors of the Hilbert space and the resulting 
DMFT equations describe a two-level system with each 
level coupled to an arbitrary external field. This problem 
(or the equivalent Anderson impurity problem) is much 
more difficult to solve than the FK impurity problem 
with a single external field coupled to an electron. 



C. Nonequilibrium Effects 

Equilibrium properties in the DMFT have been stud- 
ied extensively. Much less effort has been placed on 
nonequilibrium calculations, where the correlated many- 
body system is driven away from equilibrium by large 
amplitude external fields or driving voltages. The formal- 
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ism for D MFT in a nonequilibrium situ ation has been ad- 
dressed bv lSchmidt and Monienl l)2002f) . They considered 
the case of a Hubbard model that interacts with a large 
amplitude electromagnetic field with a constant driving 
frequency u). The formalism employs Keldysh techniques 
to derive the three Green's functions of relevance: the 
advanced, retarded, and so-called lesser Green's function. 
They find that in the steady state, one can perform cal- 
culations in frequency space, if the driving frequency is 
uniform. They applied their calculations to the Hub- 
bard model, and used the iterated perturbation theory 
to evaluate the impurity-problem dynamics. Neverthe- 
less, their formalism is general, and can be used for FK 
model problems as well. In addition, one could choose 
to perform calculations purely in the time domain. This 
would allow one to calculate the response of the system 
to large amplitude, and fast (i.e. femtosecond) electro- 
magnetic pulses. In this case, one would be interested in 
the transient response that brings the system back to the 
equilibrium state after being disturbed. 

There is another class of problems that are ideal for 
being considered with Keldysh techniques. Namely, the 
problem of the steady-state current response to a volt- 
age applied over an inhomogeneous device. One can 
imagine stacking correlated planes (described by the FK 
model) in between a semi-infinite number of noninter- 
acting metallic planes (above and below) that form bal- 
listic current leads for the device. One enforces current 
conservation from plane to plane, which then determines 
the change in the voltage from plane to plane in the de- 
vice. In this fashion, one can calculate the self-consistent 
current-voltage characteristic of a nonlinear device that 
contains correlated materials (and also the voltage pro- 
file throughout the device). The leads can be described 
by a wide class of Hamiltonians, including the Hartree- 
Fock description of a ferromagnet, or a Bardeen-Cooper- 
Schrieffer superconductor, or a diffusive FK metal, and 
the barrier material (of arbitrary width) can be described 
by the FK model tuned through the metal-insulator tran- 
sition. Then, since the total action is quadratic in the 
fermionic variables (although it is not time-translation- 
invariant) the path integral over the Keldysh contour 
can be evaluated by taking the determinant of the cor- 
responding continuous (matrix) operator. Discretizing 
the operator as was done when we described the local- 
ized electron spectral function, will then allow for the 
response to be solved exactly. We believe this will be 
the first example of a correlated electronic system that 
can undergo a metal-insulator transition, and have it's 
nonequilibrium response determined exactly. The calcu- 
lation is nontrivial, but is likely to produce a series of 
interesting results. One can apply such a formalism to a 
wide variety of devices such as Josephson junctions, spin- 
valve transistors, ballistic spin filters, or thermoelectric 
coolers. The choices are endless, as long as one always 
restricts the effective actions to remain quadratic in the 
fermionic variables. Generalization to more complicated 
systems is also possible, but usually requires one to re- 



sort to approximate methods to solve the corresponding 
impurity problems. 



VI. CONCLUSIONS 

It has been over thirty years since Falicov and Kimball 
introduced their model to describe the physics behind 
rare-earth and transition-metal metal-insulator transi- 
tions. The model has a rich history and much effort 
has been devoted to solving the model and illustrating 
its properties. The field was dramatically advanced four- 
teen years ago when Brandt and Mielsch showed that 
the model could be solved exactly in the limit of large 
dimensions. Since then much work has been done on 
investigating the solutions of the model. 

In this review we have covered the exact solution of 
the Falicov-Kimball model (and the related static Hol- 
stein model) with dynamical mean-field theory. Our fo- 
cus was on developing the formalism with a path-integral 
approach that concentrated on the Matsubara frequency 
representation of the fermions. Falicov-Kimball model 
physics is a mature field where nearly all thermody- 
namic properties of the solutions have been worked out 
in the limit of large spatial dimensions. This includes 
charge-density-wave transitions, metal-insulator transi- 
tions, phase separation, electrical and thermal transport, 
inelastic light scattering, and so on. The model has also 
been employed in more applied problems such as investi- 
gating thermodynamic properties of Josephson junctions. 

These exact solutions are useful for two reasons. First, 
there is a growing list of materials that appear to be able 
to be described by FK model physics, including YbInCu4, 
EuNi2(Sii_ 2; Ge 2 ;)2, Ta^N, and N1I2. Similar models have 
been also applied to the doped manganites and to di- 
luted magnetic semiconductors. Second, exact solutions 
are useful benchmarks for illustrating the phenomena as- 
sociated with strong correlations and for testing different 
approximation methods for their accuracy. 

The newer frontiers that will be examined in the future 
include a thorough investigation of equilibrium properties 
in finite-dimensional systems, the addition of hybridiza- 
tion or itineracy of the /-electrons and the examination 
of nonequilibrium effects in all dimensions. The success 
to date with this model, indicates that the time is ripe 
for more applications in the near future. It seems likely 
the FK model can be employed to determine a number of 
interesting applications in correlated devices, which can 
lay the groundwork for later efforts in more complicated 
correlated systems. We look forward to seeing how the 
field will develop in the ensuing years and hope that this 
review will encourage others to enter the field and rapidly 
contribute to its future development. 
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m a (rriri) z-component of spin of electrons 

2s + 1 number of a states 

t hopping 

t* rescaled hopping 

v(e) velocity operator 

w(x) phonon distribution function 

wo = 1 — wi empty site density 

wi localized electron (ion) density 

Xi phonon coordinate at site i 

r cc irreducible charge vertex 

A superconducting gap at T = 

E self energy 

£l uc unit cell volume 

a an anharmonic phonon potential (quartic) 

P 1/T 

[3 an anharmonic phonon potential (cubic) 

7(q) parameter in charge susceptibility 

?j(q) deviation from local susceptibility 

£o {£,()) parameters for F v calculation 

k spring constant 

A (2) dynamical mean field 

p chemical potential, conduction electrons 

Pb Bohr magneton 

fif chemical potential, localized electrons 

ivi — iixTH bosonic Matsubara frequency 

tp (tp) fermionic Grassman variable 

p(e) bare conduction electron DOS 

pdc dc resistivity 

p e itinerant electron density 

Pf localized electron density 

o{v) optical conductivity 

do conductivity unit 

Gdc dc conductivity 

r imaginary time 

T a (ui) relaxation time 

x' spin susceptibility 

X cc conduction electron charge susceptibility 

X cc pair-field susceptibility 

X ci mixed charge susceptibility 

X^ localized electron charge susceptibility 

Xhyb hybridization susceptibility 

Xo bare susceptibility 

Xt time-dependent conduction electron field 
iuj n = iirT(2n + 1) fermionic Matsubara frequency 

TL Hamiltonian 

C Lorenz number 

Z partition function 



